CHAPTER 


Lagrange’s Equations 


The theoretical development of the laws of motion of bodies is a problem of such interest and 
importance that it has engaged the attention of all the most eminent mathematicians since 
the invention of dynamics as a mathematical science by Galileo, and especially since the 
wonderful extension which was given to that science by Newton. Among the successors of 
those illustrious men, Lagrange has perhaps done more than any other analyst to give extent 
and harmony to such deductive researches, by showing that the most varied consequences 
respecting the motions of systems of bodies may be derived from one radical formula; the 
beauty of the methods so suiting the dignity of the results as to make of his great work a 
kind of scientific poem. 

—William Rowan Hamilton, 1834 

Armed with the ideas of the calculus of variations, we are ready to set up the version 
of mechanics published in 1788 by the Italian-French astronomer and mathematician 
Lagrange (1736-1813). The Lagrangian formulation has two important advantages 
over the earlier Newtonian formulation. First, Lagrange’s equations, unlike Newton’s, 
take the same form in any coordinate system. Second, in treating constrained systems, 
such as a bead sliding on a wire, the Lagrangian approach eliminates the forces of 
constraint (such as the normal force of the wire, which constrains the bead to remain 
on the wire). This greatly simplifies most problems, since the constraint forces are 
usually unknown, and this simplification comes at almost no cost, since we usually 
do not want to know these forces anyway. 

In Section 7.1,1 prove that Lagrange’s equations are equivalent to Newton’s second 
law for a particle moving unconstrained in three dimensions. The extension of this 
result to N unconstrained particles is surprisingly straightforward, and I leave the 
details for you to supply (Problem 7.7). In the next few sections, I take up the harder, 
and more interesting, case of constrained systems. I begin with some simple examples 
and important definitions (such as degrees of freedom). Then, in Section 7.4,1 prove 
Lagrange’s equations for a particle constrained to move on a curved surface (leaving 
the general case to Problem 7.13). Section 7.5 offers several examples, some of which 
are distinctly easier to set up in the Lagrangian formulation than in the Newtonian. In 
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Section 7.6,1 introduce the curious terminology of “ignorable coordinates.” Finally, 
after some summarizing remarks in Section 7.7, the chapter concludes with three 
sections on topics which, although very important, could be omitted on a first reading. 
In Section 7.8,1 discuss how the laws of energy and momentum conservation appear 
in Lagrangian mechanics. Section 7.9 describes how Lagrange’s equations can be 
extended to include magnetic forces, and Section 7.10 is an introduction to the idea 
of Lagrange multipliers. 

Throughout this chapter, except in Section 7.9, I treat only the case that all 
nonconstraint forces are conservative or can, at least, be derived from a potential 
energy function. This restriction can be significantly relaxed, but already includes 
most of the applications that you are likely to meet in practice. 


7.1 Lagrange’s Equations for Unconstrained Motion 


Consider a particle moving unconstrained in three dimensions, subject to a conserva¬ 
tive net force F(r). The particle’s kinetic energy is, of course, 

T = jmv 2 = |mr 2 = jm(x 2 + y 2 + z 2 ), (7.1) 

and its potential energy is 


U = U(r) = U(x 9 y 9 z). 

The Lagrangian function, or just Lagrangian, is defined as 


(7.2) 


L = T -U. 


(7.3) 


Notice first that the Lagrangian is the KE minus the PE. It is not the same as the 
total energy. You are certainly entitled to ask why the quantity T — U should be of 
any interest. There seems to be no simple answer to this question except that it is, as 
we shall see directly. Notice also that I am using a script £ for the Lagrangian 1 (to 
distinguish it from the angular momentum L and a length L) and that £ depends on the 
particle’s position (x, y, z) and its velocity (x, y, z); that is, £ = £(x, y, z, x, y, z). 

Let us consider the two derivatives, 


and 


a£ _ du 

dx 3x 


(7.4) 


3 £ 

dx 


dT 

— = mx = 
dx 


(7.5) 


1 This notation gets into difficulty in field theories where the Lagrangian is often denoted by L, 
and £ is used for the Lagrangian density , but this won’t be a problem for us. 
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Differentiating the second equation with respect to time and remembering Newton’s 
second law, F x = p x (I take for granted that our coordinate frame is inertial), we see 
that 

d A = ± d A. (7.6) 

dx dt dx 

In exactly the same way we can prove corresponding equations in y and z. Thus 
we have shown that Newton’s second law implies the three Lagrange equations (in 
Cartesian coordinates so far): 

dL d dL dL d dL J dL d dL 

— =-- , — =- r , and — =- 7 . (7.7) 

dx dt dx dy dt dy dz dt dz 

You can easily check that the argument just given works equally well in reverse, so that 
(for a single particle in Cartesian coordinates, at least) Newton’s second law is exactly 
equivalent to the three Lagrange equations (7.7). The particle’s path as determined 
by Newton’s second law is the same as the path determined by the three Lagrange 
equations. 

Our next step is to recognize that the three equations of (7.7) have exactly the 
form of the Euler-Lagrange equations (6.40). Therefore, they imply that the integral 
S = f L dt is stationary for the path followed by the particle. That this integral, called 
the action integral, is stationary for the particle’s path is called Hamilton’s principle 2 
(after its inventor, the Irish mathematician, Hamilton, 1805-1865) and can be restated 
as follows: 


H Hamilton's Principle 

The actual path which a particle follows between two points 1 and 2 in a given 
time interval, t l to f 2 , is such that the action integral 



Although we have so far proved this principle only for a single particle and in Cartesian 
coordinates, we are going to find that it is valid for a huge class of mechanical systems 
and for almost any choice of coordinates. 

So far we have proved for a single particle that the following three statements are 
exactly equivalent: 


0 . 

Try not to be confused by the unlucky circumstance that Hamilton’s principle is one possible 
statement of the Lagrangian formulation of classical mechanics (as opposed to the Hamiltonian 
formulation). 
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1. A particle’s path is determined by Newton’s second law F = ma. 

2. The path is determined by the three Lagrange equations (7.7), at least in 
Cartesian coordinates. 

3. The path is determined by Hamilton’s principle. 

Hamilton’s principle has found generalizations in many fields outside classical 
mechanics (field theories, for example) and has given a unity to various diverse areas 
of physics. In the twentieth century it has played an important role in the formulation 
of quantum theories. However, for our present purposes its great importance is that it 
lets us prove that Lagrange’s equations hold in more-or-less any coordinate system: 

Instead of the Cartesian coordinates r = (jc, y, z), suppose that we wish to use 
some other coordinates. These could be spherical polar coordinates (r, 0,0), or 
cylindrical polars (p, 0 , z), or any set of “generalized coordinates” q h q 2 , q 2 , with the 
property that each position r specifies a unique value of (q h q 2 , # 3 ) and vice versa; 
that is, 


q t = q t (r) for i = 1, 2, and 3, (7.9) 

and 

r = r(q l ,q 2 ,q 3 ). (7.10) 

These two equations guarantee that for any value of r = (x, y, z) there is a unique 
(q\,q 2 , q 3 ) and vice versa. Using (7.10) we can rewrite (x, y,z) and (x, y, z) in terms 
of ( q h q 2 , q 3 ) and (q l9 q 2 , q 3 ). Next, we can rewrite the Lagrangian £ = \mx 2 — U (r) 
in terms of these new variables as 

£ = £(#1> #3’ #2’ #3) 

and the action integral as 

rh 

S= C(q l ,q 2 ,q 3 ,q l ,q 2 ,q 3 )dt. 

Now, the value of the integral S is unaltered by this change of variables. Therefore, 
the statement that S is stationary for variations of the path around the correct path 
must still be true in our new coordinate system, and, by the results of Chapter 6, this 
means that the correct path must satisfy the three Euler-Lagrange equations, 

3£ d 3£ 3£ d 3£ , 3£ d 3£ 

— =-r, -=-and — =--, (7.11) 

dq x dt dq x dq 2 dt 3 q 2 3 q 3 dt 3 q 3 

with respect to the new coordinates q h q 2 , and q 3 . Since these new coordinates are 
any set of generalized coordinates, the qualification “in Cartesian coordinates” can be 
omitted from the statement (2) above. This result — that Lagrange’s equations have 
the same form for any choice of generalized coordinates — is one of the two main 
reasons that the Lagrangian formalism is so useful. 

There is one point about our derivation of Lagrange’s equations that is worth 
keeping at the back of your mind. A crucial step in our proof was the observation 
that (7.6) was equivalent to Newton’s second law F x = p x , which in turn is true only 
if the original frame in which we wrote down £ = T — U is inertial. Thus, although 
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Lagrange’s equations are true for any choice of generalized coordinates q h q 2 , q 2 — 
and these generalized coordinates may in fact be the coordinates of a noninertial 
reference frame — we must nevertheless be careful that, when we first write down 
the Lagrangian L = T — U, we do so in an inertial frame. 

We can easily generalize Lagrange’s equations to systems of many particles, but 
let us first look at a couple of simple examples. 



example 7.1 One Particle in Two Dimensions; Cartesian Coordinates 

Write down Lagrange’s equations in Cartesian coordinates for a particle mov¬ 
ing in a conservative force field in two dimensions and show that they imply 


Newton’s second law. (Of course, we have already proved this, but it is worth 
seeing it worked out explicitly.) 

The Lagrangian for a single particle in two dimensions is 


£ = £(jc, y,x,y) — T — U = \m{x 2 + y 2 ) - U(x, y). 
To write down the Lagrange equations we need the derivatives 

ac du t ac ar 

— =-= I\ and — = — = mx , 

dx dx dx dx 


(7.12) 

(7.13) 


with corresponding expressions for the y derivatives. Thus the two Lagrange 
equations can be rewritten as follows: 


ac 

dx 

ac 

ay 


d a£ 

«=> F x 


— 

= mx 

dt dx 



d dL 

Fy 

= my 

dt 3y 


F = ma. 




(7.14) 


Notice how in (7.13) the derivative dL/dx is the x component of the force, and 
dL/dx is the x component of the momentum (and similarly with the y components). 
When we use generalized coordinates we shall find that d&/dq i9 

although not necessarily a force component, plays a role very similar to a force. 
Similarly, dL / dq t , although not necessarily a momentum component, acts very like a 
momentum. For this reason we shall call these derivatives the generalized force and 
generalized momentum respectively; that is, 


and 


dc 


dq< 


— (ith component of generalized force) 


ac 

Hi 


= (ith component of generalized momentum). 


With these notations, each of the Lagrange equations (7.11) 


3L d 3L 


(7.15) 


(7.16) 


dq t dt dq t 



242 Chapter 7 Lagrange’s Equations 


takes the form 

(generalized force) = (rate of change of generalized momentum) (7.17) 
I shall illustrate these ideas in the next example. 

example 7.2 One Particle in Two Dimensions; Polar Coordinates 

I Find Lagrange’s equations for the same system, a particle moving in two dimen¬ 
sions, using polar coordinates. 

I 

As m all problems in Lagrangian mechanics, our first task is to write down 
the Lagrangian L = T — U in terms of the chosen coordinates. In this case we 
have been told to use polar coordinates, as sketched in Figure 7.1. This means 
| the components of the velocity are v r = r and — r0, and the kinetic energy 
| is T — ±mv 2 = \m(r 2 + r 2 0 2 ). Therefore, the Lagrangian is 

H=L(r,<p,r,<j>) = T-U=±m(r 2 + r 2 4> 2 ^ -U{r,<p). (7.18) 

Given the Lagrangian, we now have only to write down the two Lagrange 
equations, one involving derivatives with respect to r and the other derivatives 
with respect to 0. 

The r Equation 

The equation involving derivatives with respect to r (the r equation) is 

3L d d C 

dr dt dr 

;o dU d . .. 

mrcp -= — \mr)—mr . (7.19) 

dr dt 

Since —dU/dr is just F r , the radial component of F, we can rewrite the r 

I 

equation as 

j F r = m(r — r<p 2 ), (7.20) 

| 

l 

I 



i 



Figure 7.1 The velocity of a particle expressed in 
two-dimensional polar coordinates. 
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which you should recognize as F r = ma r , the r component of F = ma, first 
derived in Equation (1.48). (The term —r0 2 is the infamous centripetal accel¬ 
eration.) That is, when we use polar coordinates (r, 0), the Lagrange equa¬ 
tion corresponding to r is just the radial component of Newton’s second law. 
(Note, however, that the Lagrangian derivation avoided the tedious calculation 
of the components of the acceleration.) As we shall see directly, the 0 equation 
works a bit differently and illustrates a remarkable feature of the Lagrangian 
approach. 


The 0 Equation 

The Lagrange equation for the coordinate 0 is 


d£j _ d d£j 

30 dt 30 


or, substituting (7.18) for L, 


3 U_ 
30 


Li 2 * 

— (mr 0). 
dt 


(7.21) 


(7.22) 


To interpret this equation, we need to relate the left side to the appropriate 
component of the force F = — Vt/. This requires that we know the components 
of Vt/ in polar coordinates: 


Vt/ = 


31/. 13 U * 

—r H-0. 

dr r 30 


(7.23) 


(If you don’t remember this, see Problem 7.5.) The 0 component of the force is 

A 

just the coefficient of 0 in F = — Vt/, that is, 


r _ 131 / 

^( 1 ) A , ' 

r dtp 


Thus the left side of (7.22) is rF which is simply the torque T on the particle 
about the origin. Meanwhile, the quantity mr 2 <p on the right can be recognized 
as the angular momentum L about the origin. Therefore, the 0 equation (7.22) 
states that 


T = 


dL 
~dt ’ 


(7.24) | 


the familiar condition from elementary mechanics, that torque equals the rate 
of change of angular momentum. 


The result (7.24) illustrates a wonderful feature of Lagrange’s equations, that when 
we choose an appropriate set of generalized coordinates the corresponding Lagrange 
equations automatically appear in a corresponding, natural form. When we choose 
r and 0 for our coordinates, the 0 equation turns out to be the equation for angular 
momentum. In fact, the situation is even better than this. Recall that I introduced the 
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notion of generalized force and generalized momentum in (7.15) and (7.16). In the 
present case, the 0 component, of the generalized force is just the torque, 

3 £ 

(0 component of generalized force) = — = T (torque) (7.25) 

30 

and the corresponding component of the generalized momentum is 

3£ 

(0 component of generalized momentum) = —- = L (angular momentum). (7.26) 

d(j) 

With the “natural” choice for the coordinates (r and 0) the 0 components of the 
generalized force and momentum turn out to be the corresponding “natural” quantities, 
the torque and the angular momentum. 

Notice that the generalized “force” does not necessarily have the dimensions of 
force, nor the generalized “momentum” those of momentum. In the present case, 
the generalized force (0 component) is a torque (that is, force x distance) and the 
generalized momentum is an angular momentum (momentum x distance). 

This example illustrates another feature of Lagrange’s equations: The 0 component 
3L/30 of the generalized force turned out to be the torque on the particle. If the 
torque happens to be zero, then the corresponding generalized momentum 3L/30 
(the angular momentum, in this case) is conserved. Clearly this is a general result: 
The ith component of the generalized force is dL/dq^. If this happens to be zero, then 
the Lagrange equation 

3£ __ d_ 3L 
3 q t dt 3 q t 

says simply that the ith component 3£/3 q t of the generalized momentum is constant. 
That is, if £ is independent of q h the ith component of the generalized force is zero, 
and the corresponding component of the generalized momentum is conserved. In 
practice, it is often easy to spot that a Lagrangian is independent of a coordinate q h 
and, if you can, then you immediately know a corresponding conservation law. We 
shall return to this point in Section 7.8. 


Several Unconstrained Particles 

The extension of the above ideas to a system of N unconstrained particles (a gas of 
N molecules, for instance) is very straightforward, and I shall leave you to fill in the 
details (Problems 7.6 and 7.7). Here I shall just sketch the argument for the case of 
two particles, mainly to show the form of Lagrange’s equations for N > 1. For two 
particles, the Lagrangian is defined (exactly as before) as £ = T — U, but this now 
means that 


£(r b r 2 , r h r 2 ) = \m^ + \m 2 x\ - U (r b r 2 ). (7.27) 

As usual, the forces on the two particles are F x = — V X U and F 2 = —V 2 U. Newton’s 
second law can be applied to each particle and yields six equations, 

F\x Plxi F\y P\y-> * ’ ^ 2 z Plz' 
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Exactly as in Equation (7.7), each of these six equations is equivalent to a correspond¬ 
ing Lagrange equation 


3 £ _ d 3£ dL_d_ ac 

dxi dt dx x ’ dyi dt dyi ’ 


3 £ d 3 £ 

3 z 2 dt 3z 2 


(7.28) 


These six equations imply that the integral S = £ dt is stationary. Finally, we can 

change to any other suitable set of six coordinates q h q 2 , • - , q^- The statement that S 
is stationary must also be true in this new coordinate system, and this implies in turn 
that Lagrange’s equations must be true with respect to the new coordinates: 

3£ _ d 3£ 3£_d_3£ 3£ _ d_ 3£ 

3#! <7/ 3^ 3^2 <7^ 3# 2 3g 6 dt dq 6 

An example of a set of six such generalized coordinates that we shall use repeatedly 
in Chapter 8 is this: In place of the six coordinates of r x and r 2 , we could use the 
three coordinates of the CM position R = ( m l r l + m 2 r 2 )/(m 1 + m 2 ) and the three 
coordinates of the relative position r = r { — r 2 . We shall find that this choice of 
coordinates leads to a dramatic simplification. For now, the main point is simply that 
Lagrange’s equations are automatically true in their standard form (7.29) with respect 
to these new, generalized coordinates. 

The extension of these ideas to the case of N unconstrained particles is entirely 
straightforward, and I leave it for you to check. (See Problem 7.7.) The upshot is that 
there are 3 N Lagrange equations 


3 £ d 3 £ 

dq t dt dcfa 


[i = 1, 2,..., 3N], 


valid for any choice of the 3N coordinates q h • • •, q 3N needed to describe the N 
particles. 


7.2 Constrained Systems; an Example 


Perhaps the greatest advantage of the Lagrangian approach is that it can handle systems 
that are constrained so that they cannot move arbitrarily in the space that they occupy. 
A familiar example of a constrained system is the bead which is threaded on a wire — 
the bead can move along the wire, but not anywhere else. Another example of a very 
constrained system is a rigid body, whose individual atoms can only move in such a 
way that the distance between any two atoms is fixed. Before I discuss the nature of 
constraints in general, I shall discuss another simple example, the plane pendulum. 

Consider the simple pendulum shown in Figure 7.2. A bob of mass m is fixed to a 
massless rod, which is pivoted at O and free to swing without friction in the vy plane. 
The bob moves in both the x and y directions, but it is constrained by the rod so 
that y/x 2 + y 2 = l remains constant. In an obvious sense, only one of the coordinates 
is independent (as x changes, the variation of y is predetermined by the constraint 
equation), and we say that the system has only one degree of freedom. One way to 
express this would be to eliminate one of the coordinates, for instance by writing 
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o 


Figure 7.2 A simple pendulum. The bob of mass m is con¬ 
strained by the rod to remain at distance / from O . 


y — V/ 2 — x 2 and expressing everything in terms of the one coordinate x. Although 
this is a perfectly legitimate way to proceed, a simpler way is to express both x and y in 
terms of the single parameter 0, the angle between the pendulum and its equilibrium 
position, as shown in Figure 7.2. 

We can express all the quantities of interest in terms of 0. The kinetic energy is 
T — \mv 2 = \ml 2 <j) 2 . The potential energy is U = mgh where h denotes the height of 
the bob above its equilibrium position and is (as you should check) h =1(1 — cos 0). 
Thus the potential energy is U = mgl( 1 — cos </>), and the Lagrangian is 

£ = T - U = \ml 2 4> 2 - mgl{\ - cos<j>). (7.30) 


Whichever way we choose to proceed — to write everything in terms of x (or y) or 
<fi — the Lagrangian is expressed in terms of a single generalized coordinate q and its 
time derivative q , in the form L = L(q, q). Now, it is a fact (which I shall not prove 
just yet) that once the Lagrangian is written in terms of this one variable (for a system 
with one degree of freedom), the evolution of the system again satisfies Lagrange’s 
equation (just as we proved for an unconstrained particle in the previous section.) 
That is, 


d£j d d£j 
3 q dt 3 q 


(7.31) 


If we choose the angle (j) as our generalized coordinate, then Lagrange’s equation 
reads 


3L _ d 3L 
3 (p dt 3 <p 


(7.32) 


The Lagrangian L is given by (7.30), and the needed derivatives are easily evaluated 
to give 


—mgl sin 0 


— (ml 2 <p) = ml 2 (j). 
dt 


(7.33) 
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Referring to Figure 7.2 you can see that the left side of this equation is just the torque 
T exerted by gravity on the pendulum, while the term ml 2 is the pendulum’s moment 
of inertia /. Since 0 is the angular acceleration a, we see that Lagrange’s equation for 
the simple pendulum simply reproduces the familiar result T = la. 


7.3 Constrained Systems in General _ 

Generalized Coordinates 

Consider now an arbitrary system of N particles, a = 1, • • •, N with positions r^. 
We say that the parameters q u • • •, q n are a set of generalized coordinates for the 
system if each position r^ can be expressed as a function of q h • • •, q n , and possibly 
the time t, 


r a = r a(q 1> -"iqmt) [« = 1, , 2V], (7.34) 

and conversely each can be expressed in terms of the r^ and possibly t , 

q t = q t ( r l9 , r N ,t) [i = 1, • • •, n\ (7.35) 

In addition, we require that the number of the generalized coordinates (n) is the 
smallest number that allows the system to be parametrized in this way. In our three- 
dimensional world, the number n of generalized coordinates for N particles is cer¬ 
tainly no more than 3 N and, for a constrained system, is usually less — sometimes 
dramatically so. For example, for a rigid body, the number of particles N may be of 
order 10 23 , whereas the number of generalized cordinates n is 6 (three coordinates to 
give the position of the center of mass and three to give the orientation of the body). 

To illustrate the relation (7.34), consider again the simple pendulum of Figure 7.2. 
There is one particle (the bob) and two Cartesian coordinates (since the pendulum is 
restricted to two dimensions). As we saw, there is just one generalized coordinate, 
which we took to be the angle 0. The analog of (7.34) is 

r = (x, y) = (/ sin0, l cos 0) (7.36) 

and expresses the two Cartesian coordinates x and y in terms of the one generalized 
coordinate 0. 

The double pendulum shown in Figure 7.3 has two bobs, both confined to a plane, 
so it has four Cartesian coordinates, all of which can be expressed in terms of the two 
generalized coordinates (j) l and 0 2 . Specifically, if we put our origin at the suspension 
point of the top pendulum, 

I-J = (Zj sin 0 l9 l x cos 0j) = r 1 (0 1 ) (7.37) 


and 


r 2 = (/i sin 0j + Z 2 sin0 2 ,Z 1 cos0 1 + / 2 cos0 2 ) = r 2 (0 1? 0 2 ). (7.38) 


Notice that the components of r 2 depend on 0j and 0 2 . 
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Figure 7.3 The positions of both masses in a double pendulum 
are uniquely specified by the two generalized coordinates 0! and 
0 2 , which can themselves be varied independently. 


In these two examples, the transformation between the Cartesian and the general¬ 
ized coordinates did not depend on the time t , but it is easy to think of examples in 
which it does. Consider the railroad car shown in Figure 7.4, which has a pendulum 
suspended from its ceiling and is being forced 3 to accelerate with a fixed accelera¬ 
tion a. It is natural to specify the position of the pendulum by the angle 0 as usual, 
but we must recognize that, in the first instance, this gives the pendulum’s position 
relative to the accelerating, and hence non-inertial, reference frame of the car. If we 
wish to specify the bob’s position relative to an inertial frame, we can choose a frame 
fixed relative to the ground, and we can easily express the position relative to this 
inertial frame in terms of the angle 0. The position of the point of suspension relative 
to the ground is (if we choose our axes and origin properly) just x s = \at 2 , and the 
position of the bob is then easily seen to be 

r = (x, y ) = (/ sin0 + ^ at 2 , / cos0) = r(0, t). (7.39) 



Figure .7.4 A pendulum is suspended from the roof of a rail¬ 
road car that is being forced to accelerate with a fixed, known 
acceleration a. 


3 The word “forced” is often used to describe a motion that is imposed by some outside agent and 
is unaffected by the internal motions of the system. In the present example, the “forced” acceleration 
of the car is assumed to be the same whatever the oscillations of the pendulum. 
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The relation between r and the generalized coordinate 0 depends on the time t , a 
possibility that I allowed for when writing (7.34). 

We shall sometimes describe a set of coordinates q h • • • ,q n as natural if the 
relation (7.34) between the Cartesian coordinates r a and the generalized coordinates 
does not involve the time t. We shall find certain convenient properties of natural 
coordinates that do not generally apply to coordinates for which (7.34) does involve 
the time. Fortunately, as the name implies, there are many problems for which the 
most convenient choice of coordinates is also natural. 4 


Degrees of Freedom 

The number of degrees of freedom of a system is the number of coordinates that 
can be independently varied in a small displacement — the number of independent 
“directions” in which the system can move from any given initial configuration. For 
example, the simple pendulum of Figure 7.2 has just one degree of freedom, while the 
double pendulum of Figure 7.3 has two. A particle that is free to move anywhere in 
three dimensions has three degrees of freedom, while a gas comprised of N particles 
has 3 N. 

When the number of degrees of freedom of an A-particle system in three dimen¬ 
sions is less than 3A, we say that the system is constrained. (In two dimensions, the 
corresponding number is 2A of course.) The bob of a simple pendulum, with one 
degree of freedom, is constrained. The two masses of a double pendulum, with two 
degrees of freedom, are constrained. The A atoms of a rigid body have just six degrees 
of freedom and are certainly constrained. Other examples are a bead constrained to 
slide on a fixed wire and a particle constrained to move on a fixed surface in three 
dimensions. 

In all of the examples I have given so far, the number of degrees of freedom 
was equal to the number of generalized coordinates needed to describe the system’s 
configuration. (The double pendulum has two degrees of freedom and needs two 
generalized coordinates, and so on.) A system with this natural-seeming property is 
said to be holonomic. 5 That is, a holonomic system has n degrees of freedom and 
can be described by n generalized coordinates, q h • • •, q n . Holonomic systems are 
easier to treat than nonholonomic, and in this book I shall restrict myself to holonomic 
systems. 

You might imagine that all systems would be holonomic, or at least that nonholo¬ 
nomic systems would be rare and bizarrely complicated. In fact, there are some quite 
simple examples of nonholonomic systems. Consider, for instance, a hard rubber ball 
that is free to roll (but not to slide nor to spin about a vertical axis) on a horizontal 
table. Starting at any position (.x, y) it can move in only two independent directions. 
Therefore, the ball has two degrees of freedom, and you might well imagine that its 


4 Natural coordinates are sometimes called scleronomous, and those that are not natural, 
rheonomous. I shall not use these outstandingly forgettable names. Nonnatural coordinates are 
also sometimes called forced, since a time dependence in the relation (7.34) is usually associated 
with a forced motion, such as the forced acceleration of the car in Figure 7.4. 

5 Many different definitions of “holonomic” can be found, not all of which are exactly equivalent. 



Chapter 7 Lagrange’s Equations 



Figure 7.5 The right triangle OPQ lies in the xy plane with sides 
OP and PQ of length c. If you roll a ball of circumference c around 
OPQ , it will return to its starting point with a changed orientation. 


configuration could be uniquely specified by two coordinates, x and y, of its center. 
But consider the following: Let us place the ball at the origin O and make a mark on 
its top point. Now, carry out the following three moves. (See Figure 7.5.) Roll the ball 
along the x axis for a distance equal to the circumference c , to a point P, where the 
mark will once again be on the top. Now roll it the same distance c in the y direction to 
Q , where the mark is again on top. Finally roll it straight back to the origin along the 
hypotenuse of the triangle OPQ. Since this last move has length V2c, it brings the ball 
back to its starting point, but with the mark no longer on the top. The position (x, y) 
has returned to its initial value, but the ball now has a different orientation. Evidently 
the two coordinates (x, y) are not enough to specify a unique configuration. In fact, 
three more numbers are needed to specify the orientation of the ball, and we need five 
coordinates in all to specify the configuration completely. The ball has two degrees 
of freedom but needs five generalized coordinates. Evidently it is a nonholonomic 
system. 

Although nonholonomic systems certainly exist, they are more complicated to an¬ 
alyze than holonomic systems, and I shall not discuss them further. For any holonomic 
system with generalized coordinates qi,---,q n and potential energy U (q\, • • • ,q n ,t) 
(which may depend on the time t as described in Section 4.5), the evolution in time 
is determined by the n Lagrange equations 


3£ _ d_ 3£ 
3 q t dt dcfa 


n], 


(7.40) 


where the Lagrangian £ is defined as usual to be £ = T — U. I shall prove this result 
in Section 7.4. 


7.4 Proof of Lagrange’s Equations with Constraints 


We are now ready to prove Lagrange’s equations for any holonomic system. To keep 
things reasonably simple, I shall treat explicitly the case that there is just one particle. 
(The generalization to arbitrary numbers of particles is fairly straightforward — see 
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Problem 7.13.) To be definite, I shall suppose the particle is constrained to move on a 
surface. 6 This means that it has two degrees of freedom and can be described by two 
generalized coordinates q x and q 2 that can vary independently. 

We must recognize that there are two kinds of forces on the particle (or particles, in 
the general case). First, there are the forces of constraint: For a bead on a wire the con¬ 
straining force is the normal force of the wire on the bead; for our particle, constrained 
to move on a surface, it is the normal force of the surface. For the atoms in a rigid 
body, the constraining forces are the interatomic forces that hold the atoms in place 
within the body. In general, the forces of constraint are not necessarily conservative, 
but this doesn’t matter. One of the objectives of the Lagrangian approach is to find 
equations that do not involve the constraining forces, which we usually don’t want to 
know anyway. (Notice, however, that if the constraining forces are nonconservative, 
Lagrange’s equations in the simple unconstrained form of Section 7.1 certainly do not 
apply.) I shall denote the net constraining force on the particle by F cstr , which in our 
case is just the normal force of the surface to which the particle is confined. 

Second, there are all the other “nonconstraint” forces on the particle, such as 
gravity. These are the forces with which we are usually concerned in practice, and 
I shall denote their resultant by F. I shall assume that the nonconstraint forces all 
satisfy at least the second condition for conservatism, so that they are derivable from 
a potential energy, U (r, t), and 


F = — VC/(r, t). (7.41) 

(If all the nonconstraint forces are actually conservative, then U is independent of t, 
but we don’t need to assume this.) The total force on our particle is F tot = F cstr + F. 
Finally, I shall define the Lagrangian, as usual, to be 

L = T - U. (7.42) 

Since U is the potential energy for the nonconstraint forces only, this definition of £ 
excludes the constraint forces. This correctly reflects that Lagrange’s equations for a 
constrained system cleverly eliminate the constraint forces, as we shall see. 


The Action Integral is Stationary at the Right Path 

Consider any two points iq and r 2 through which the particle passes at times t x and 
t 2 . I shall denote by r (t) the “right” path, the actual path that our particle follows 
between the two points, and by R (t) any neighboring “wrong” path between the same 
two points. It is convenient to write 

R(0 - r(f) + 6(0, (7.43) 


6 Actually, it is a bit hard to imagine how to constrain a particle to a single surface so that it 
can’t jump off. If this worries you, you can imagine the particle sandwiched between two parallel 
surfaces with just enough gap between them to let it slide freely. 
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which defines e(t) as the infinitesimal vector pointing from r(r) on the right path to 
the corresponding point R(0 on the wrong path. Since I shall assume that both of 
the points r(7) and R (t) lie in the surface to which the particle is confined, the vector 
e(t) is contained in the same surface. Since both r(r) and R(7) go through the same 
endpoints, e(t) = 0 at t x and t 2 - 

Let us denote by 5 the action integral 

rh 

S= / £(R,R,0<*, (7.44) 

Jt x 

taken along any path R (t) lying in the constraining surface, and by S 0 the correspond¬ 
ing integral taken along the right path r(Y). As I shall now prove, the integral S is 
stationary for variations of the path R(0, when R(f) = r(0 or, equivalently, when 
the difference e is zero. Another way to say this is that the difference in the action 
integrals 


8S = S - S 0 (7.45) 

is zero to first order in the distance € between the paths, and this is what I shall prove. 

The difference (7.45) is the integral of the difference between the Lagrangians on 
the two paths, 


8L = £(R, R, t) - £(r, r, t). (7.46) 

If we substitute R (t) = r (t) + e(t) and 

L(r, r, t) = T — U — \mr 2 — U (r, t), 

this becomes 7 

= \m [V + kj 2 - r 2 j -[U(r + €,t)-U (r, 0] 

= mr4 - (-VL + 0(6 2 ), 

where O (e 2 ) denotes terms involving squares and higher powers of 6 and k. Returning 
to the difference (7.45) in the two action integrals, we find that, to first order in 6, 

f 2 8£j dt = f 2 [mr • k — € -VU]dt. (7.47) 

Jt\ Jt x 

The first term in the final integral can be integrated by parts. (Recall that this just 
means moving the time derivative from one factor to the other and changing the sign.) 
The difference e is zero at the two endpoints, so the endpoint contribution is zero, and 
we get 

8S = - [ 2 €-[mr + VU]dt. (7.48) 

Jh 


7 To understand the second term in the second line, recall that /(r + e) — /(r) ^ € • V/, for 
any scalar function /(r). See Section 4.3. 
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Now, the path r (t) is the “right” path and satisfies Newton’s second law. Therefore the 
term mr is just the total force on the particle, F tot = F cstr + F. Meanwhile Vt/ = —F. 
Therefore, the second term in (7.48) cancels the second piece of the first, and we are 
left with 

SS = - f 2 €-F cstT dt . (7.49) 

Jt\ 

But the constraint force F cstr is normal to the surface in which our particle moves, 
while e lies in the surface. Therefore e • F cstr = 0, and we have proved that 8S ~ 0. 
That is, the action integral is stationary at the right path, as claimed. 8 


The Final Proof 

We have proved Hamilton’s principle, that the action integral is stationary at the path 
which the particle actually follows. However, we have proved it, not for arbitrary 
variations of the path, but rather for those variations of path that are consistent with 
the constraints — that is, paths which lie in the surface to which our particle is 
constrained. This means that we cannot prove Lagrange’s equations with respect to 
the three Cartesian coordinates. On the other hand, we can prove them with respect to 
the appropriate generalized coordinates. We are assuming that our particle is confined 
by holonomic constraints to move on a surface, that is, a two-dimensional subset 
of the full three-dimensional world. This means that the particle has two degrees of 
freedom and can be described by two generalized coordinates, q x and q 2 , that can be 
varied independently. Any variation of q x and q 2 is consistent with the constraints. 9 
Accordingly, we can rewrite the action integral in terms of q x and q 2 as 

rh 

s = I L(q h q 2 ,q x ,q 2 ,t)dt, (7.50) 

A 

and this integral is stationary for any variations of q x and q 2 about the correct path 
lq x (t ), q 2 {t) ]. Therefore, by the argument of Chapter 6 the correct path must satisfy 
the two Lagrange equations 

dc ' d dL J ac d ac 

— =-- and — =-- . (7.51) 

dq x dt dq x dq 2 dt dq 2 

The proof that I have given here applies directly only to a single particle in three 
dimensions, constrained to move on a two-dimensional surface, but the main ideas 
of the general case are all present. The generalization is, for the most part, relatively 


8 The observation that the integrand in (7.49) is zero is really the crucial step in our proof. When 
you consider the generalization of the proof to an arbitrary constrained system (for instance, if you 
do Problem 7.13), you will find that there is a corresponding step and that the corresponding term 
is zero, for the same reason: The forces of constraint would do no work in a displacement that is 
consistent with the constraints. Indeed this is one possible definition of a force of constraint. 

9 For example, if our surface is a sphere, centered at the origin, then the generalized coordinates 
q i, q 2 could be the two angles 0, <j> of spherical polar coordinates. Any variation of 0 and <fi is 
consistent with the constraint that the particle remain on the sphere. 
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straightforward (see Problem 7.13), and meanwhile I hope that I have said enough 
to convince you of the truth of the general result: For any holonomic system, with n 
degrees of freedom and n generalized coordinates, and with the nonconstraint forces 
derivable from a potential energy U(q h • • •, q n , t ), the path followed by the system 
is determined by the n Lagrange equations 


i 


r* i 

«M«w* | ; 
*'***** J. 4 


(7.52) 


where L is the Lagrangian L = T — U and U = U(q x , • • •, q n , t) is the total potential 
energy corresponding to all the forces excluding the forces of constraint. 

It was essential to our proof of Lagrange’s equations that the nonconstraint forces 
be conservative (or, at a minimum, that they satisfy the second condition for conser¬ 
vatism) so that they are derivable from a potential energy, F = — Vf/. If this is not 
true, then Lagrange’s equations may not hold, at least in the form (7.52). An obvi¬ 
ous example of a force that does not satisfy this condition is sliding friction. Sliding 
friction cannot be regarded as a force of constraint (it is not normal to the surface) 
and cannot be derived from a potential energy. Thus, when sliding friction is present, 
Lagrange’s equations do not hold in the form (7.52). Lagrange’s equations can be 
modified to include forces like friction (see Problem 7.12), but the result is clumsy 
and I shall confine myself to situations where the equations (7.52) do hold. 


7.5 Examples of Lagrange’s Equations 


In this section I present five examples of the use of Lagrange’s equations. The first two 
are sufficiently simple that they can be easily solved within the Newtonian formalism. 
My main purpose for including them is just to give you experience with using the 
Lagrangian approach. Nevertheless, even these simple cases show some advantages 
of the Lagrangian over the Newtonian formalism; in particular, we shall see how the 
Lagrangian approach obviates any need to consider the forces of constraint. The last 
three examples are sufficiently complex that solution using the Newtonian approach 
requires considerable ingenuity; by contrast, the Lagrangian approach lets us write 
down the equations of motion almost without thinking. 

The examples given here illustrate an important point to recognize about La¬ 
grange’s equations: The Lagrangian formalism always (or nearly always) gives a 
straightforward means of writing down the equations of motion. On the other hand, it 
cannot guarantee that the resulting equations are easy to solve. If we are very lucky, 
the equations of motion may have an analytic solution, but, even when they do not, 
they are the essential first step to understanding the solutions and they often suggest 
a starting point for an approximate solution. The equations of motion can give simple 
answers to certain subsidiary questions. (For instance, once we have the equations of 
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motion, we can usually find the positions of equilibrium of a system very easily.) And 
we can always solve the equations of motion numerically for given initial conditions. 

The Lagrangian method is so important that it certainly deserves more than just 
five examples. However, the crucial thing is that you work through several examples 
yourself; therefore I have given plenty of problems at the end of the chapter, and it is 
essential that you work several of these as soon as possible after reading this section. 


example 7.3 Atwood’s Machine 

Consider the Atwood machine first met in Figure 4.15 and shown again in Figure 
7.6, in which the two masses m x and m 2 are suspended by an inextensible string j 
(length /) which passes over a massless pulley with frictionless bearings and 
radius R. Write down the Lagrangian £, using the distance x as generalized co¬ 
ordinate, find the Lagrange equation of motion, and solve it for the acceleration 
x. Compare your results with the Newtonian solution. j 

Because the string has fixed length, the heights x and y of the two masses j 
cannot vary independently. Rather, x + y + ttR = /, the length of the string, so | 
that y can be expressed in terms of x as j 

y = —x + const. (7.53) j 

Therefore, we can use x as our one generalized coordinate. From (7.53) we see | 
that y — —x, so that the kinetic energy of the system is 

T — ^m x x 2 + \m 2 y 2 — \(yn x + m 2 )x 2 , | 



Figure 7.6 An Atwood machine consisting of two masses, m x 
and ra 2 , suspended by a massless inextensible string that passes 
over a massless, frictionless pulley of radius R. Because the 
string’s length is fixed, the position of the whole system can 
be specified by a single variable, which we can take to be the 
distance x. 
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I while the potential energy is 

I 

j JJ — —m \gx — m 2 gy = —(mi — m 2 )gx + const. 


Combining these, we find the Lagrangian 

& = T — U = ^ (m! + m 2 )x 2 + (m 1 — m 2 )gx , (7.54) 

where I have dropped an uninteresting constant. 

The Lagrange equation of motion is just 

d£j d d £j 
dx dt dx 

or, substituting (7.54) for £, 

{m x -m 2 )g - {m x +m 2 )x, (7.55) 

which we can solve at once to give the desired acceleration 


x = 


m ] — m 2 
m \ + m 2 


(7.56) 


By choosing m x and m 2 fairly close together, one can make this acceleration 
much less than g, and hence much easier to measure. Therefore, the Atwood 
machine gave an early and reasonably accurate method for measuring g. 

The corresponding Newtonian solution requires us to write down Newton’s 
second law for each of the masses separately. The net force on m x is m x g — F t 
where F t is the tension in the string. (This is the force of constraint and needed 
no consideration in the Lagrangian solution.) Thus Newton’s second law for m x 
is 


rriig — F t = mpc. 


In the same way, Newton’s second law for m 2 reads 
J F t - m 2 g = m 2 x. 

I 

| (Remember that the upward acceleration of m 2 is the same as the downward 
acceleration of m j.) We see that the Newtonian approach has given us two equa¬ 
tions for two unknowns, the required acceleration x and the force of constraint 
F t . By adding these two equations, we can eliminate F t and arrive at precisely 
the equation (7.55) of the Lagrangian method and thence the same value (7.56) 
for x. 

1 

The Newtonian solution of the Atwood machine is too simple for us to get 
very excited by an alternative solution. Nevertheless, this simple example does 
illustrate how the Lagrangian approach allows us to ignore the unknown (and 
usually uninteresting) force of constraint and to eliminate at least one step of 
j the Newtonian solution. 
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example 7.4 A Particle Confined to Move on a Cylinder 

Consider a particle of mass m constrained to move on a frictionless cylinder of 
radius R , given by the equation p = R in cylindrical polar coordinates (p , 0, z), 
as shown in Figure 7.7. Besides the force of constraint (the normal force of the 
cylinder), the only force on the mass is a force F = —kr directed toward the | 
origin. (This is the three dimensional version of the Hooke’s-law force.) Using | 
z and 0 as generalized coordinates, find the Lagrangian £. Write down and solve j 
Lagrange’s equations and describe the motion. 

Since the particle’s coordinate p is fixed at p = R, we can specify its position | 
by giving just z and 0, and since these two coordinates can vary independently 
the system has two degrees of freedom and we can use (z, 0) as generalized 
coordinates. The velocity has v p = 0, = R<p , and v z = z. Therefore the kinetic | 

energy is 

T = ^mv 2 = ^m(R 2 <p 2 + z 2 ). 

The potential energy for the force F = —kr is (Problem 7.25) U = \kr 2 , where r 
is the distance from the origin to the particle, given by r 2 = R 2 + z 2 (see Figure 
7.7). Therefore 


t/ = i*(J? 2 + z 2 ), | 

and the Lagrangian is 

£ = \m{R 2 0 2 + i 2 ) - ifc(i? 2 + z 2 ). (7.57) 1 

| 

Since the system has two degrees of freedom, there are two equations of j 
motion. The z equation is j 



Figure 7.7 A mass m is confined to the surface of the cylinder 
p = R and subject to a Hooke’s law force F = — kr. 
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The 0 equation is even simpler. Since £ does not depend on 0, it follows that 
3£/30 — 0 and the 0 equation is just 


3£ d 3£ 
30 dt 30 


r/ 9 • 

or 0 = — mRd). 
dt 


(7.59) 


The z equation (7.58) tells us that the mass executes simple harmonic motion in 
the z direction, with z = A cos (cot — 8). The 0 equation (7.59) tells us that the 
quantity mR 2 <p is constant, that is, that the angular momentum about the z axis 
is conserved — a result we could have anticipated since there is no torque in 
this direction. Because p is fixed, this implies simply that 0 is constant, and the 
mass moves around the cylinder with constant angular velocity 0, at the same 
time that it moves up and down in the z direction in simple harmonic motion. 


These two examples illustrate the steps to be followed in solving any problem by 
the Lagrangian method (provided all constraints are holonomic and the nonconstraint 
forces are derivable from a potential energy, as we are assuming): 

1. Write down the kinetic and potential energies and hence the Lagrangian 
£ = T — U , using any convenient inertial reference frame. 

2. Choose a convenient set of n generalized coordinates q h • • •, q n and find 
expressions for the original coordinates of step 1 in terms of your chosen 
generalized coordinates. (Steps 1 and 2 can be done in either order.) 

3. Rewrite £ in terms of q l9 • • •, q n and q l9 • • •, q n . 

4. Write down the n Lagrange equations (7.52). 

As we shall see, these four steps provide an almost infallible route to the equations of 
motion of any system, no matter how complex. Whether the resulting equations can 
be easily solved is another matter, but even when they cannot, just having them is a 
huge step toward understanding a system and an essential step to finding approximate 
or numerical solutions. 

The next two examples illustrate how the Lagrangian approach can give the equa¬ 
tions of motion, almost effortlessly, for problems that would require considerable care 
and ingenuity using Newtonian methods. 


example 7.5 A Block Sliding on a Wedge 

Consider the block and wedge shown in Figure 7.8. The block (mass m) is free 
to slide on the wedge, and the wedge (mass M) can slide on the horizontal table, 
both with negligible friction. The block is released from the top of the wedge, 

j with both initially at rest. If the wedge has angle a and the length of its sloping 

| 

face is Z, how long does the block take to reach the bottom? 

The system has two degrees of freedom, and a good choice of the two 
generalized coordinates is, as shown, the distance q x of the block from the top 
of the wedge and the distance q 2 of the wedge from any convenient fixed point 
on the table. The quantity we need to find is the acceleration q x of the block 

relative to the wedge, since with this we can quickly find the time required to 

i 

! slide the length of the wedge. Our first task is to write down the Lagrangian, and 
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Figure 7.8 A block of mass m slides down a wedge of mass M , 
which is free to slide over the horizontal table. 

it is often safest to do this in Cartesian coordinates, and then rewrite it in terms 
of the chosen generalized coordinates. 

The kinetic energy of the wedge is just T M = \Mq 2 2 but that of the block 
is more complicated. The block’s velocity relative to the wedge is q x down the 
slope, but the wedge itself has a horizontal velocity q 2 relative to the table. The 
velocity of the block relative to the inertial frame of the table is the vector sum of 
these two. Resolving into rectangular components (x to the right, y downward), 
we find for the velocity of the block relative to the table 

v = (v x , v y ) = (cfa cos ot + q 2 , q x sin o'). 

Thus the kinetic energy of the block is 

10 0 1 0 0 
T m — ^ m(v x + v) = 2 m (q l + q 2 + 2q x q 2 cos a). 

(I used the identity cos 2 a + sin 2 a = 1 to simplify this.) The total kinetic energy 
of the system is 

T = T m + T m = \{M + m)q^ + \m(qf + 2q x q 2 cosa). (7.60) 

The potential energy of the wedge is a constant, which we may as well take to 
be zero. That of the block is — mgy , where y = q x sin a is the height of the block 
measured down from the top of the wedge. Therefore 

U = —mgq x sin a 


and the Lagrangian is 

1 0 1 • 0 

= T — U = 2 (M + m)g 2 + + 2 q x q 2 cos a) + mgq x sin o'. (7.61) 

Once we have found the Lagrangian in terms of the generalized coordinates 
q x and q 2 , all we have to do is to write down the two Lagrange equations, one 
for q x and one for q 2 , and then solve them. The q 2 equation (which is a little 
simpler) is 

9L d 3L 


dq 2 dt dq 2 


(7.62) 
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but, since L in (7.61) is clearly independent of q 2 , this just tells us that the 
generalized momentum dL/dq 2 is constant, 


Mq 2 + m(q 2 + q x cos a) = const 


(7.63) 


— a result you will recognize as conservation of the total momentum in the x 
direction (and something you could have written down without any help from 
Lagrange). 

The q x equation 


3L d d£j 
dq x dt dq x 


(7.64) 


is more complicated, since neither derivative is zero. Sustituting (7.61) for £, 
we can write this as 


mg sin a — —m (^1 + ^2 cos a ) 
dt 


m(q x + q 2 cos a) 


(7.65) 


Differentiating (7.63) we see that 


q 2 = - —- q ! cos o', 

M + m 

which lets us eliminate q 2 from (7.65) and solve for q x \ 


(7.66) 



• • 

<h 


g sin o' 
m cos 2 a 
M + m 


(7.67) 


Armed with this value for q x we can quickly answer the original question: Since 
the acceleration down the slope is constant, the distance traveled down the slope 
in time t is \q x t 2 , and the time to travel the length / is just y/2l/q\, with q x given 
by (7.67). More interesting than this answer is to check that the formula (7.67) for 
q x agrees with common sense in various special cases. For example, if a = 90°, 
(7.67) implies that q x = g , which is clearly right; and, if M -> 00 , (7.67) implies 
that q x —> g sin a, which is the well-known acceleration for a block on a fixed 
incline and clearly makes sense. I leave it as an exercise (Problem 7.19) to 
check that in the limit that M —> 0, our answers agree with what you could 
have predicted. 


example 7.6 A Bead on a Spinning Wire Hoop 

A bead of mass m is threaded on a frictionless circular wire hoop of radius R. 
) The hoop lies in a vertical plane, which is forced to rotate about the hoop’s 
vertical diameter with constant angular velocity <p = &>, as shown in Figure 
7.9. The bead’s position on the hoop is specified by the angle 6 measured up 

| 

I from the vertical. Write down the Lagrangian for the system in terms of the 









Section 7.5 Examples of Lagrange’s Equations 


261 



Figure 7.9 A bead is free to move around the frictionless wire 
hoop, which is spinning at a fixed rate co about its vertical axis. 
The bead’s position is specified by the angle #; its distance from 
the axis of rotation is p = R sin #. 


generalized coordinate # and find the equation of motion for the bead. Find any 
equilibrium positions at which the bead can remain with # constant, and explain 
their locations in terms of statics and the “centrifugal force” mco 2 p (where p 
is the bead’s distance from the axis). Use the equation of motion to discuss the 
stability of the equilibrium positions. 

Our first task is to write down the Lagrangian. Relative to a nonrotating 
frame, the bead has velocity RO tangential to the hoop and pco — (R sin 0)(o 
normal to the hoop (the latter due to the spinning of the hoop with angular 


velocity co). Thus the kinetic energy is T — \mv l = \mR z (6 A + or sin #). 
The gravitational potential energy is easily seen to be U — mgR( 1 — cos#), 
measured from the bottom of the hoop. Therefore, the Lagrangian is 


L = jmR 2 (9 2 + co 2 sin 2 #) — mgR{ 1 — cos#), 


(7.68) 


and the Lagrange equation is 


d dfl 


or mR 2 co 2 sin# cos# — mgR sin# = mR L 9 


2 n 


a# dt a# 


Dividing through by mR 2 , we arrive at the desired equation of motion: 


6 = (nr cos# — g/R) sin#. 


(7.69) 


Although this equation cannot be solved analytically in terms of elementary 
functions, it can, nevertheless, tell us lots about the system’s behavior. To 
illustrate this, let us use (7.69) to find the equilibrium positions of the bead. 
An equilibrium point is any value of # — call it # 0 — satisfying the following 
condition: If the bead is placed at rest (# = 0) at # = # 0 , then it will remain at 
rest at # 0 . This condition is guaranteed if 9 = 0. (To see this, note that if 9 = 0, 
then 9 doesn’t change and remains zero, which means that # doesn’t change 
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and remains equal to 0 o ) Thus to find the equilibrium positions we have only to 
equate the right side of (7.69) to zero: 


( co 2 cos 9 - g/R) sin 9 = 0. (7.70) 

1 

This equation is satisfied if either of the two factors is zero. The factor sin 0 is 
zero if 6 = 0 or n . Thus the bead can remain at rest at the bottom or top of the 
hoop. The first factor in (7.70) vanishes when 


cos 0 = 


8 

co 2 R 


| Since | cos#| must be less than or equal to 1, the first factor can vanish only 

9 

when co > g/R. When this condition is satisfied, there are two more equilibrium 
positions at 


0 o = ± arccos 



(7.71) 


We conclude that when the hoop is rotating slowly (co 2 < g/R), there are just 
two equilibrium positions, at the bottom and top of the hoop, but when it rotates 
fast enough (co 2 > g/R), there are two more, symmetrically placed on either 
side of the bottom, as given by (7.71). 10 

Perhaps the simplest way to understand the various equilibrium positions is 
in terms of the “centrifugal force.” In most introductory physics courses, the 
centrifugal force is dismissed as an abomination to be avoided by all right- 
thinking physicists. As long as we confine our attention to inertial frames, this 
is a correct (and certainly a safe) point of view. Nevertheless, as we shall see 
in Chapter 9, from the point of view of a noninertial rotating frame there is a 
perfectly real centrifugal force mco 2 p (perhaps more familiar as mv 2 /p ), where 
p is the object’s distance from the axis of rotation. Thus, taking the point of 
view of a fly perched on the rotating hoop, we can understand the equilibrium 
positions as follows: At the bottom or top of the hoop, the bead is on the axis of 
rotation and p — 0; therefore, the centrifugal force mco 2 p is zero. Furthermore, 
the force of gravity is normal to the hoop, so there is no force tending to move 
the bead along the wire and the bead remains at rest. The other two equilibrium 
points are a little subtler: At any position off the axis (such as that shown in 
Figure 7.9) the centrifugal force is nonzero and has a component pushing the 

| 

bead outward along the wire; meanwhile the force of gravity has a component 

| 

pulling the bead inward along the wire (provided the bead is below the halfway 
marks, 9 — ±7t/2). Ateitherofthepoints given by (7.71), these two components 
are balanced (check this for yourself— Problem 7.28) and the bead can remain 

I 

i at rest. 

An equilibrium point 0 o is not especially interesting unless it is stable — that 
is, the bead, if nudged a little away from 0 o , moves back toward 6 0 . Using our 


10 Notice that when co 2 = g/R the two extra positions given by (7.71) have just come into 
existence and coincide with the first point at the bottom with 6 — ±0. 
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equation of motion (7.69), we can easily address this issue, and I’ll start with 
the equilibrium at the bottom, 0 — 0. As long as # remains close to 0, we can 
set cos # ~ 1 and sin 0 & # and approximate the equation by 

6 = {co 2 — g/R)9 [6» near 0]. (7.72) j 

If the hoop is rotating slowly {co 1 < g/R), this has the form 

6 = (negative number)#. 

If we nudge the bead to the right (0 > 0), then since # is positive # is negative, 
and the bead accelerates to the left, that is, back toward the bottom. If we nudge 
it to the left, (# < 0), then # becomes positive, and the bead accelerates to the j 
right, which is again back toward the bottom. Either way, the bead returns toward 
the equilibrium, which is, therefore, stable. 

If we speed up the rotation of the hoop, so that co 2 > g/R, then the approxi¬ 
mate equation of motion (7.72) takes the form | 

0 = (positive number)#. | 

Now a small displacement to the right makes # positive, and the bead accelerates | 
away from the bottom. Similarly a displacement to the left makes 6 negative, j 
and again the bead accelerates away from the bottom. Thus, as we increase co | 
past the critical value where co 2 = g/R, the equilibrium at the bottom changes 
from stable to unstable. j 

The equilibrium at the top (# = n) is alway unstable (see Problem 7.28). 
This is easy to understand from our discussion of the centrifugal force. Near the j 
top of the hoop, both the centrifugal and gravitational forces tend to push the 
bead away from the top, so there is no chance of a restoring force to pull it back I 
to the equilibrium position. 

The other two equilibrium positions only exist when co 2 > g/R , and are 
easily seen to be stable: The equation of motion (7.69) is 

I 

<9 = (cn 2 cos <9 - g/f?)sin6>. (7.73) | 

To be definite, let us consider the equilibrium on the right with 0 < # < jt/2. At 
the equilibrium point, the term in parenthesis on the right of (7.73) is zero, while 
sin # is positive. If we increase # a little (bead moves up and to the right), sin # 
remains positive, but the term in parenthesis becomes negative. (Remember, 
cos # is a decreasing function in this quadrant.) Thus # becomes negative, and 
the bead accelerates back toward its equilibrium point. If we decrease # a little | 
from the equilibrium, then # becomes positive, and again the bead accelerates 
back toward equilibrium. Therefore the equilibrium on the right is stable. As you 
would expect, a similar analysis shows that the same is true of the equilibrium 
on the left. | 

We arrive at the following interesting story: When the hoop is rotating slowly 
(co 2 < g/R), there is just one stable equilibrium, at # = 0. If we speed up 
the rotation, then as co passes the critical value where co 2 = g/R, this original 
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equilibrium becomes unstable, but two new stable equilibrium points appear, 
emerging from 6=0 and moving out to the right and left as we increase co 
still more. This phenomenon — the disappearance of one stable equilibrium and 
the simultaneous appearance of two others diverging from the same point — is 

| 

| called a bifurcation and will be one of our principal topics in Chapter 12 on 
chaos theory. 

It is interesting to note that the device of this example was used by James 
Watt (1736-1829) as a governor for his steam engines. The device rotated with 
the engine, and as the engine sped up the bead rose on the hoop. When the 
angular velocity co reached some predetermined maximum allowable value, the 
bead, arriving at a corresponding height, caused the supply of steam to be shut 



This example illustrates another strength of the Lagrangian method that was 
mentioned back in Section 7.1: The generalized coordinates can even be coordinates 
relative to a noninertial reference frame, just as long as the frame in which the 
Lagrangian L = T — U was orginally written was inertial. In this example, the angle 
6 was the polar angle of the bead, measured in the noninertial rotating frame of the 
hoop, but the Lagrangian (7.68) was defined as L = T — U with T and U evaluated 
in the inertial frame relative to which the hoop rotates. 11 

In the next and final example of this section, we pursue the previous example of 
the bead on the rotating hoop, and obtain approximate solutions of the equation of 
motion in the neighborhood of the stable equilibrium points. 


example 7.7 Oscillations of the Bead near Equilibrium 

Consider again the bead of the previous example and use the equation of motion 
to find the bead’s approximate behavior in the neighborhood of the stable 
equilibrium positions. 

When co 2 < g/R, the only stable equilibrium is at the bottom of the hoop 
with 6 = 0. As long as 6 remains small, we can approximate the equation of 
motion (7.73) by setting sin 6 ^ 6 and cos 6 ^ 1 to give 

6 = —(g/R — co 2 )0 [6 near 0] 

I = —Q 2 6 (7.74) 

where the second line introduces the frequency 


| Q = yJg/R - CO 2 . 

I 

I 

As long as co 2 < g/R, this defines Q as a real positive number, and we recog¬ 
nize (7.74) as the equation for simple harmonic motion with frequency £2. We 


11 Example 7.5 was another instance: The coordinate q x gave the position of the block relative 
to the accelerating frame of the wedge, but the kinetic energy T was evaluated in the inertial frame 
of the table. For another example, see Problem 7.30. 
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conclude that a bead which is displaced a little from the stable equilibrium at | 
0=0, executes harmonic motion with frequency Q, 

0(t) = A cos(Qt — 8). (7.75) j 

If we speed up the rate of the hoop’s rotation until co 2 > g/R, then Q becomes 
pure imaginary, and, since cos ia = cosh a, our solution (7.75) becomes a hy¬ 
perbolic cosine, which grows with time, correctly reflecting that the equilibrium 
at 0 = 0 has become unstable. 

Once co 2 > g/R , there are two stable equilibrium positions given by (7.71) 
and located symmetrically to the right and left of the bottom of the hoop. As 
you might expect, these behave in the same way, and to be definite I shall focus ] 
on the one on the right. Let us denote its position by 0 = 6 0 , where, according 
to (7.70), 0 o satisfies | 


co 2 cos6 0 — g/R = 0. (7.76) 

Let us now imagine the bead placed close to 0 o at 

0 = 0 o + € 

and investigate the time dependence of the small parameter e. Once again we 
can approximate the equation of motion (7.73), though this requires more care. 
If we approximate the factors of cos(# 0 + e) and sin(0 o + e) by the first two 
terms of their Taylor series, 


cos(0 o + e) ^ cos 0 o — e sin 0 o and sin(0 o + c) ^ sin 0 o + c cos 0 o (7.77) 


then the equation of motion (7.73) becomes 

0 = [co 2 cos(# 0 + c) — g/R] sin(0 o + e) [0 near 0 O ] 

= [co 2 cos 0 o — 6 co 2 sin 6 0 — g/R ] [sin 0 o + 6 cos 0 o ]. (7.78) 

By (7.76) the first and third terms in the first square bracket cancel, leaving just 
the middle term —eco 2 sin 0 O . To lowest order in e we can drop the second term 
of the second bracket, and, since 0 is the same as €, we are left with 


€ = -€ co 2 sin 2 0 o = -Q ,2 €. (7.79) 

Here the second equality defines the frequency = co sin 6 0 , or, using (7.76), 


Q f = 



(7.80) 


(see Problem 7.26). Equation (7.79) is the equation for simple harmonic motion. 
Therefore, the parameter e oscillates about zero, and the bead itself oscillates 
about the equilibrium position 0 o with frequency Q f . 
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7.6 Generalized Momenta and Ignorable Coordinates 


As I have already mentioned, for any system with n generalized coordinates q t 
(i — n), we refer to the n quantities 3£/3 q t — F t as generalized forces and 

3 £/dqi = p t as generalized momenta. With this terminology, the Lagrange equation, 


can be rewritten as 


3£ d 3£ 

3 q t dt dcfo ’ 




(7.81) 


(7.82) 


That is, “generalized force = rate of change of generalized momentum.” In particular, 
if the Lagrangian is independent of a particular coordinate q h then F ( = 3£ /dq { = 0 
and the corresponding generalized momentum p t is constant. 

Consider, for example, a single projectile subject only to the force of gravity. 
The potential energy is U — mgz (if we use Cartesian coordinates with z measured 
vertically up), and the Lagrangian is 


£ = £(x, y, z, x, y, z) = \m(x 2 + y 2 + z 2 ) — mgz. (7.83) 

With respect to Cartesian coordinates, the generalized force is just the usual force 
(3£ /dx = —dU/dx = F x , etc.) and the generalized momentum is just the usual mo¬ 
mentum (3£/3i = mx = p x , etc.) Because £ is independent of x and y, it immedi¬ 
ately follows that the components p x and p y are constant, as we already knew. 

In general, the generalized forces and momenta are not the same as the usual 
forces and momenta. For instance, we saw in Equations (7.25) and (7.26) that in two- 
dimensional polar coordinates the (j> component of the generalized force is the torque, 
and that of the generalized momentum is actually the angular momentum. In any case, 
when the Lagrangian is independent of a coordinate q t the corresponding generalized 
momentum is conserved. Thus, if the Lagrangian of a two-dimensional particle is 
independent of </>, then the particle’s angular momentum is conserved — another im¬ 
portant result (and one that is clear from the Newtonian perspective as well). When 
the Lagrangian is independent of a coordinate q t , that coordinate is sometimes said 
to be ignorable or cyclic. Obviously it is a good idea, when possible, to choose co¬ 
ordinates so that as many as possible are ignorable and their corresponding momenta 
are constant. In fact, this is perhaps the main criterion in choosing generalized coor¬ 
dinates for any given problem: Try to find coordinates as many as possible of which 
are ignorable. 

We can rephrase the result of the last three paragraphs by noting that the statement 
“£ is independent of a coordinate qf is equivalent to saying “£ is unchanged, or 
invariant , when q t varies (with all the other # • held fixed).” Thus we can say that if 
£ is invariant under variations of a coordinate q t then the corresponding generalized 
momentum p t is conserved. This connection between invariance of £ and certain 
conservation laws is the first of several similar results relating invariance under 
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transformations (translations, rotations, and so on) to conservation laws. These results 
are known collectively as Noether’s theorem, after the German mathematician Emmy 
Noether (1882-1935). I shall return to this important theorem in Section 7.8. 


7.7 Conclusion 


The Lagrangian version of classical mechanics has the two great advantages that, 
unlike the Newtonian version, it works equally well in all coordinate systems and 
it can handle constrained systems easily, avoiding any need to discuss the forces of 
constraint. If the system is constrained, one must choose a suitable set of independent 
generalized coordinates. Whether or not there are constraints, the next task is to write 
down the Lagrangian £ in terms of the chosen coordinates. The equations of motion 
then follow automatically in the standard form 


3£ _ d_dL_ 
dq t dt dq t 


[/ = !,•••, n\. 


There is, of course, no guarantee that the resulting equations will be easy to solve, 
and in most real problems they are not, requiring numerical solution or at least some 
approximations before they can be solved analytically. 

Even in problems that are only moderately complicated, like the examples of 
Section 7.5, finding the equations of motion by Lagrange’s method is remarkably 
easier than by using Newton’s second law. Indeed, some purists object that the 
Lagrangian approach makes life too easy, removing the need to think about the 
physics. 

The Lagrangian formalism can be extended to include more general systems than 
those considered so far. One important case is that of magnetic forces, which I take up 
in Section 7.9. Dissipative forces, such as friction or air resistance, can sometimes be 
included, but it should be admitted that the Lagrangian formalism is primarily suited 
to problems where dissipative forces are absent or, at least, negligible. 

The final three sections of this chapter treat three advanced topics, all of which 
are centrally important in Lagrangian mechanics, but all of which could be omitted 
on a first reading. In Section 7.8, I give two more examples of the remarkable 
connection between invariance under certain transformations and conservation laws. 
This connection, known as Noether’s theorem, is important in all of modern physics, 
but especially in quantum physics. Section 7.9 discusses how to include magnetic 
forces in Lagrangian mechanics, another topic of great importance in quantum theory. 
Finally, Section 7.10 introduces the method of Lagrange multipliers. This technique 
appears in many different guises in many areas of physics, but I shall restrict myself to 
some simple examples in Lagrangian mechanics. These last three sections are arranged 
to be self-contained and independent. You could study all of them, none of them, or 
any selection in between. 
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7.8 More about Conservation Laws* 


* The material of this section is more advanced than the preceding sections, and you should 
feel free to omit it on a first reading. Be aware, however, that the material discussed here is 
needed before you read Section 11.5 and Chapter 13. 

In this section I shall discuss how the laws of conservation of momentum and energy 
fit into the Lagrangian formulation of mechanics. Since we derived the Lagrangian 
formulation from the Newtonian, anything that we already knew about conservation 
laws, based on Newtonian mechanics, will naturally still be true in Lagrangian me¬ 
chanics. Nevertheless, we can gain some new insights by examining the conservation 
laws from a Lagrangian perspective. Furthermore, much modem work takes the La¬ 
grangian formulation (based on Hamilton’s principle, for example) as its starting point. 
In this context, it is important to know what can be said about conservation laws strictly 
within the Lagrangian framework. 


Conservation of Total Momentum 

We already know from Newtonian mechanics that the total momentum of an isolated 
system of N particles is conserved, but let us examine this important property from the 
Lagrangian point of view. One of the most prominent features of an isolated system 
is that it is translationally invariant; that is, if we transport all N particles bodily 
through the same displacement e, nothing physically significant about the system 
should change. This is illustrated in Figure 7.10, where we see that the effect of moving 
the whole system through the fixed displacement e is to replace every position by 


r l rj + €, r 2 -> r 2 + *, • • •, r N -*r N + €. (7.84) 

In particular, the potential energy must be unaffected by this displacement, so that 
U(r l + €, • • •, r N + €, t) = U (r l5 • • •, r N ,t) (7.85) 



Figure 7.10 An isolated system of N particles is translationally 
invariant, which means that when every particle is transported 
through the same displacement e, nothing physically significant 
changes. 
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or, more briefly, 


8U = 0 


where 8U denotes the change in U under the translation (7.84). Clearly the velocities 
are unchanged by the translation (7.84). (Adding a constant € to all the doesn’t 
change the r a .) Therefore 8T = 0, and hence 


8L=0 


(7.86) 


under the translation (7.84). This result is true for any displacement e. If we choose € 
to be an infinitesimal displacement in the x direction, then all of the x coordinates 
x 1? • • •, x N increase by 6, while the y and z coordinates are unchanged. For this 
translation, the change in £ is 


xr _ dL f 

8£j — 6 — + 


3x 


, 3£ A 

_J_ £ - - 0. 


3x 


N 


This implies that 


N 


T — = o. 

* J a x 

a=l 01 

Now using Lagrange’s equations we can rewrite each derivative as 


(7.87) 


9L d dL 


d 


dx 


a 


dt 3x^ dt 


Pax 


'a 


where p ax is the x component of the momentum of particle a. Thus (7.87) becomes 


d_ _ d_ p _ n 

^dt Pax ~dt x ~ 

a=\ 

where P x is the x component of the total momentum P = p a . By choosing the 
small displacement € successively in the y and z directions, we can prove the same 
result for the y and z components, and we reach the conclusion that, provided the 
Lagrangian is unchanged by the translation (7.84), the total momentum of the N- 
particle system is conserved. This connection between translational invariance of £ 
and conservation of total momentum is another example of Noether’s theorem. 


Conservation of Energy 

Finally, I would like to discuss conservation of energy from the Lagrangian point of 
view. The analysis proves somewhat complicated, but introduces a number of ideas 
that are important in more advanced work, particularly in the Hamiltonian formulation 
of mechanics (Chapter 13). 

As time advances the function L(q h • • •, q n > 4b ' * * > 4m 0 changes, both because 
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t is changing, and because the q J s and g’s change with the evolving system. Thus, by 
the chain rule, 


d . . 3 Z . a 3C .. 3 Z 

Xj(^i, * • •, q n , cfa, * ? 0 / ~ Qi T~ / . cfo 

“ 9^- V a?; 




9* 


(7.88) 


Now, by Lagrange’s equation, I can replace the derivative in the first sum on the right 
by 


9L 


Hi 


d dZ 
dt 3 q t 


d 




Meanwhile, the derivative in the second sum on the right of (7.88) is just the general¬ 
ized momentum p t . Thus, I can rewrite (7.88) as 


d £ = J2 (pm + p &)+ dL 


dt 


3 1 


d \—\ / . \ dZ 

= jX m + 


dt 


(7.89) 


Now, for many interesting systems, the Lagrangian does not depend explicitly on 
time; that is, dZ/dt = 0. When this is the case, the second term on the right of (7.89) 
vanishes. If we move the left side of (7.89) over to the right, we see that the time 
derivative of the quantity p^cp — Z is zero. This quantity is so important that it has 
its own symbol, 


, K = Y J Pi<li-t, (7.90) 

/ = 1 

and is called the Hamiltonian of the system. With this terminology, we can state the 
following important conclusion: 


If the Lagrangian Z does not depend explicitly on time (that is, dZ/dt = 0), 
then the Hamiltonian jC is conserved. 


The discovery of any conservation law is a momentous event and is enough to justify 
saying that the Hamiltonian is an important quantity. In fact, it goes much further than 
this. As we shall see in Chapter 13, the Hamiltonian 3i is the basis of the Hamiltonian 
formulation of mechanics, in just the same way that Z is the basis of Lagrangian 
mechanics. 

For the moment, the chief importance of our newly discovered Hamiltonian is that 
in many situations it is in fact just the total energy of the system. Specifically, we shall 
prove that, provided the relation between the generalized coordinates and Cartesians 
is time-independent , 


r a **a(*7l’ * * ’ ’ 9Vi)’ 


(7.91) 



Section 7.8 More about Conservation Laws* 


271 


the Hamiltonian J~C is just the total energy, 

3T = T + U. (7.92) 

You may recall that we agreed in Section 7.3 to describe generalized coordinates 
that satisfy (7.91) as natural ; thus, we can paraphrase the result (7.92) to say that, 
provided the generalized coordinates are natural, 3i is just the total energy T + U. 
To prove this, let us express the total kinetic energy T = \ m a*a * n terms °f the 
generalized coordinates q h ■ • •, q n . First, differentiating (7.91) with respect to t and 
using the chain rule, we find that 12 


. V- dr a . 

r « = l^ * 

7^ a ?.- 


(7.93) 


If we now form the scalar product of this equation with itself, we find 


r 2 = y 

a / v 

j 

where I have renamed the summation indices as j and k to avoid future confusion. 
The kinetic energy is now given as a triple sum, which I can reorganize and write as 13 



T = m m o?l = 2 A jk qjq k (7.94) 

a j,k 

where Aj k is shorthand for the sum 


Ajk — Aj k (q i, • • •» 



(7.95) 


We can now evaluate the generalized momentum p { by differentiating (7.94) with 
respect to q t (Problem 7.45), 


Pi = 



(7.96) 


Returning to Equation (7.90) for the Hamiltonian, we can rewrite the sum on the 
right as 


(7.97) 


If the relation (7.91) were explicitly time-dependent there would be one extra term in this 
expression for r a , namely dr a /dt. This extra term would invalidate the conclusion (7.98) below that 
< K = T + U. 

1 T 

We can restate the result (7.94) to say that, provided the generalized coordinates are natural, 
the kinetic energy T is a homogeneous quadratic function of the generalized velocities q t . This result 
plays an important role in several later developements. See, for instance, Section 11.5. 


pkn = = J2 A iMi = 2T 

i i V j ' i,j 
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where the last step follows from (7.94). Therefore 

% = Piki - £ = 2T - (T - U) = T + U. (7.98) 

i 

That is, provided the transformation between the Cartesian and generalized coordi¬ 
nates is time-independent, as in (7.91), the Hamiltonian J~C is just the total energy of 
the system. 

I have already proved that, provided the Lagrangian is independent of time, the 
Hamiltonian is conserved. Thus we now see that time independence of the Lagrangian 
[together with the condition (7.91)] implies conservation of energy. We can rephrase 
the time independence of £ by saying that £ is unchanged by translations of time, 
*—►> + €. Thus the result we have just proved is that invariance of £ under time 
translations is related to energy conservation, in much the same way that invariance 
of £ under translations of space (r -+ r -b e) is related to conservation of momentum. 
Both results are manifestations of Noether’s famous theorem. 


7.9 Lagrange’s Equations for Magnetic Forces* 


* This section requires a knowledge of the scalar and vector potentials of electromagnetism. 
Although the ideas described here play an important role in the quantum-mechanical treatment 
of magnetic fields, they will not be used again in this book. 

Although I have so far consistently defined the Lagrangian as L = T — U, there are 
systems, such as a charged particle in a magnetic field, which can be treated by the 
Lagrangian method, but for which £ is not just T — U. The natural question to ask 
is then: What is the definition of the Lagrangian for such systems? This is the first 
question I address. 


Definition and Nonuniqueness of the Lagrangian 

Probably the most satisfactory general definition of a Lagrangian for a mechanical 
system is this: 



In other words, a Lagrangian is any function £ for which Lagrange’s equations are 
true for the system under consideration. 
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Obviously, for the systems that we have discussed so far, the old definition £ = 
T — U fits this new definition. But the new definition is much more general. In 
particular, it is easy to see that our new definition does not define a unique Lagrangian 
function. For example, consider a single particle in one dimension and suppose that 
we have found a Lagrangian £ for this particle. That is, the equation of motion of the 
particle is 


ac _ d_ ac 

dx dt dx 


(7.99) 


Now let f(x,x) be any function for which 


dx dt dx 


(7.100) 


(It is easy to think up such a function, for instance, / = xx.) If we replace £ in 
(7.99) by 

£' = £ + / 


then, by virtue of (7.100), the Lagrangian £' gives exactly the same equation of motion 
as £. 

The lack of uniqueness of the Lagrangian is similar to, though more radical than, 
the familiar lack of uniqueness in the potential energy (to which one can add any 
constant without changing any of the physical predictions). The crucial point is that 
any function £ which gives the right equation of motion has all of the features that we 
require of a Lagrangian (for instance, that the integral / £ dt is stationary at the right 
path) and so is just as acceptable as any other such function £. If, for a given system, 
we can spot a function £ that leads to the right equation of motion, then we don’t 
need to debate whether it is the “right” Lagrangian — if it gives the right equation of 
motion, then it is just as right as any other conceivable Lagrangian. 


Lagrangian for a Charge in a Magnetic Field 

Consider now a particle (mass m and charge q) moving in electric and magnetic fields 
E and B. The force on the particle is the well-known Lorentz force F = q (E + v x B), 
so Newton’s second law reads 

rar = g(E + r x B). (7.101) 

To reformulate (7.101) in Lagrangian form, we have only to spot a function £ for 
which the three Lagrange equations are the same as (7.101). This can be done using 
the scalar and vector potentials, V (r, t) and A(r, t), in terms of which the two fields 
can be written 14 


E = -VV - 


dX 

~dt 


and B = V x A. 


(7.102) 


14 See, for example, David J. Griffiths, Introduction to Electrodynamics , (Prentice-Hall, 1999), 
p. 416-417. 
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I now claim that the Lagrangian function 15 

L(r, r, t) = \mr 2 — q(V — r • A) (7.103) 

= jm(x 2 + y 2 + z 2 ) — q(V — xA x — yA y — zA z ) (7.104) 

has the desired property, that it reproduces Newton’s second law (7.101). To check 
this, let us examine the first of the three Lagrange equations, 


3L d 3L 
dx dt dx 


(7.105) 


To see what this implies we have to evaluate the two derivatives of the proposed 
Lagrangian (7.104): 


ac 

dx 


= -q 



.3 A. 
— x—- 
dx 


.3 Ay .3 A z \ 

y —- - z—- 
dx dx ) 


(7.106) 


and 

3L 

—- = mx + qA x . 
dx 

When we differentiate this with respect to t, we must remember that A x = 
A x (x , y, z, t). As t varies, x, y, z move with the particle and, by the chain rule, we 
find 


_r/_3L 
dt dx 


= mx + q 



+ y- 


.3 A. 


dy 


. dA 

+ z-- + 

3 z 



(7.107) 


Substituting (7.106) and (7.107) into (7.105), cancelling the two terms in i, and 
rearranging, we find that Lagrange’s equation (the x component, with the proposed 
Lagrangian) is the same as 


mx — —q 


dV dA. 

+ 


dx 


dt 


+ qy 


3 A 3 
dx 


or, according to (7.102), 




(7.108) 

dx J 


mx — q (E x + yB z - zB y ) 


(7.109) 


which you will recognize as the x component of Newton’s second law (7.101). Since 
the y and z components work in the same way, we conclude that Lagrange’s equations, 
with the proposed Lagrangian (7.104), are exactly equivalent to Newton’s second law 
for the charged particle. That is, we have successfully recast Newton’s second law for 
a charged particle into Lagrangian form with the Lagrangian (7.104). 

Using the Lagrangian (7.104), one can solve various problems involving charged 
particles in electric and magnetic fields. (See Problem 7.49 for an example.) Theoret- 


15 Notice that you can, if you wish, write this as L = T — U, but U is certainly not the usual 
PE, since it depends on the velocity r; U is sometimes called a “velocity-dependent PE,” but notice 
that it is not true that the force on the charge is —VC/. 
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ically, the most important conclusion of this analysis emerges when we evaluate the 
generalized momentum. For example, 

ai: • , . 

Px = — = mx + qA x . 

OX 

Since the y and z components work the same way, we conclude that 

(generalized momentum, p) = m\ + qA. (7.110) 

That is, the generalized momentum is the mechanical momentum m\ plus a magnetic 
term qA. This result is at the heart of the quantum theory of a charged particle in 
a magnetic field, where it turns out that the generalized momentum corresponds to 
the differential operator —ifiN (where h is Planck’s constant over 2tt), so that the 
quantum analog of the mechanical momentum m\ is the operator —ihV — qA. 


7.10 Lagrange Multipliers and Constraint Forces* 


* The method of Lagrange multipliers is used in many areas of physics. Nevertheless, we shan y t 
be using it again in this book, and you could omit this section without any loss of continuity. 

In this section, we discuss the method of Lagrange multipliers. This powerful method 
finds application in several areas of physics and takes on quite different appearances 
in different contexts. Here I shall treat only its application to Lagrangian mechanics, 16 
and to keep the analysis simple I shall restrict the discussion to two-dimensional 
systems with just one degree of freedom. 

We have seen that one of the strengths of Lagrangian mechanics is that it can bypass 
all of the forces of constraint. However, there are situations where one actually needs 
to know these forces. For example, the designer of a roller coaster needs to know the 
normal force of the track on the car to know how strong to build the track. In this 
case, we can still use a modified form of Lagrange’s equations, but the procedure is 
somewhat different: We do not choose generalized coordinates q h • • •, q n all of which 
can be independently varied. (Remember that it was the independence of q h • • •, q n 
that let us use the standard Lagrange equations without worrying about constraints.) 
Instead we use a larger number of coordinates and use Lagrange multipliers to handle 
the constraints. 

To illustrate this procedure, we’ll consider a system with just two rectangular 
coordinates x and y, which are restricted by a constraint equation of the form 17 

fix, y) = const. (7.111) 

For example, we could consider a simple pendulum with just one degree of freedom 
(as in Figure 7.2). In treating this by the standard Lagrange approach one would use 


16 For applications to other kinds of problems, see, for example, Mathematical Methods in the 
Physical Sciences by Mary Boas (Wiley, 1983), Ch. 4, Section 9 and Ch. 9, Section 6. 

\T 

We’ll see directly by example that some typical constraints can be put in this form. In fact, it 
is fairly easy to show that any holonomic constraints can be. 
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the one generalized coordinate 0, the angle between the pendulum and the vertical, 
and avoid any discussion of the constraints. If, instead, we choose to use the original 
rectangular coordinates x and y, then we must recognize that these coordinates are 
not independent; they satisfy the constraint equation 


f(x,y) = Vx 2 + y 2 = i 


where l is the length of the pendulum. We shall find that the method of Lagrange 
multipliers lets us accomodate this constraint, determine the time dependence of x and 
y, and find the tension in the rod. As a second example, consider the Atwood machine 
of Figure 7.6. In our previous treatment we used the one generalized coordinate x, 
the position of the mass m h but we could instead use both coordinates x and y (the 
positions of both masses), provided we remember that the constancy of the string’s 
length imposes the constraint that 

f(x, y) — x + y = const. 

Here too, Lagrange multipliers will let us accomodate this constraint, solve for the 
time dependence of x and y, and find the constraint force, which is here the tension 
in the string. 

To set up our new method we start from Hamilton’s principle. Our Lagrangian has 
the form £ (x, x, y, y). (We could allow it to have an explicit dependence on t as well, 
but to simplify notation I shall assume it doesn’t.) The proof of Hamilton’s principle 
given in Section 7.4 applies even when the coordinates are constrained. (Indeed it 
was designed to allow for constraints.) Thus we can conclude as before that the action 
integral 


2 

S= I £(x, x, y, y) dt (7.112) 

Jt\ 

is stationary when taken along the actual path followed. If we denote this “right” path 
by x(0, y(t) and imagine a small displacement to a neighboring “wrong” path, 


xt ;\ " x( !\ + (’-ns) 

y(t) -» y(t) + 8y(t ) 

then, provided the displacement is consistent with the constraint equation, the action 
integral (7.112) is unchanged, 8S = 0. To exploit this, we must write 8S in terms of 
the small displacements 8x and 8y : 

8S= [ (—8x + —8x + —8y + —8y) dt. (7.114) 

J \dx 3x dy dy J 

The second and fourth terms can be integrated by parts, and we conclude that 


8S = 



9£ d 3£ 

-) Sx dt + 

3x dt dx 



3 £ d 3£ 


3y dt dy 


8ydt = 0 (7.115) 


for any displacements 8x and Sy consistent with the constraints. 
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If (7.115) were true for any displacements, we could immediately prove two 
separate Lagrange equations, one for jc and one for y. (Choosing 8y = 0, we would be 
left with just the first integral; since this has to vanish for any choice of Sx, the factor 
in parentheses would have to be zero, which implies the usual Lagrange equation with 
respect to jc. And similarly for y.) This is exactly the correct conclusion for the case 
that there are no constraints. 

However, there are constraints, and (7.115) is only true for displacements Sx and 
Sy consistent with the constraints. Therefore, we proceed as follows: Since all points 
with which we are concerned satisfy /(jc, y) = const, the displacement (7.113) leaves 
f(x,y) unchanged, so 


8f = -±-8x + —8y - 0 (7.116) 

dx 3 y 

for any displacement consistent with the constraints. Since this is zero, we can multiply 
it by an arbitrary function X(t) —this is the Lagrange multiplier — and add it to 
the integrand in (7.115), without changing the value of the integral (namely zero). 
Therefore 


35- 


+ 


/( 

/( 


dL w 3/ d 3L . , 

-h X(t) -I 8x dt 

dx dx dt dx 

3L ,3/ d 3£t e J A 

- h k(t) -I 8y dt — 0 

dy 3 y dt 3 y 


(7.117) 


for any displacement consistent with the constraints. Now comes the supreme cunning: 
So far X(t) is an arbitrary function of t, but we can choose it so that the coefficient of 
8x in the first integral is zero. That is, by choice of the multiplier A,(r), we can arrange 
that 



along the actual path of the system. This is the first of two modified Lagrange equations 
and differs from the usual equation only by the extra term involving A on the left. With 
the multiplier chosen in this way, the whole first integral in (7.117) is zero. Therefore 
the second integral is also zero (since their sum is), and this is true for any choice 
of 8y. (The constraint places no restriction on 8x or 8y separately — it only fixes 8x 
once 8y is chosen, or vice versa.) Therefore the coefficient of 8y in this second integral 
must also be zero and we have a second modified Lagrange equation, 



( 7 * 119 ) 
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We now have two modified Lagrange equations for the two unknown functions 
x(t) and y(t). This elegant result has been bought at the price of introducing a third 
unknown function, the Lagrange multiplier X(t). To find three unknown functions, we 
need three equations, but fortunately a third equation is already at hand, the constraint 
equation 


/(*, y) 355 const 


The three equations (7.118), (7.119), and (7.120) are sufficient, in principle at least, to 
determine the coordinates x(t) and y(t) and the multiplier X(t). Before we illustrate 
this with an example, there is one more bit of theory to develop. 

So far the Lagrange multiplier X(t) is just a mathematical artifact, introduced to 
help us solve our problem. However, it turns out to be closely related to the forces of 
constraint. To see this, we have only to look more closely at the modified Lagrange 
equations (7.118) and (7.119). The Lagrangian of our present discussion has the 
form 


L = \m x x 2 + \m 2 y 2 - U(x, y ). 


(In a problem like the simple pendulum, x and y are the two coordinates of a single 
mass and m l = m 2 . In a problem like the Atwood machine, there are two separate 
masses and m x and m 2 are not necessarily equal.) Inserting this Lagrangian into 
(7.118), we find 


at/ 

dx 


, , 3 / 

+ a— = m,x. 
dx 


(7.121) 


Now, on the left side —dU/dx is the x component of the nonconstraint force. (Re¬ 
member that U was defined as the potential energy of the nonconstraint forces.) 
On the right, m l x is the x component of the total force, equal to the sum of 
the nonconstraint and constraint forces. Thus m{x — —dU/dx + F^ str . Canceling 
the term —dU/dx from both sides of (7.121), we reach the important conclusion 
that 


k— = F cstr (7.122) 

dx x 

with a corresponding result for the y components. This then is the significance of the 
Lagrange multiplier: Multiplied by the appropriate partial derivatives of the constraint 
function /(x, y), the Lagrange multiplier X(t) gives the corresponding components 
of the constraint force. 

Let us now see how these ideas work in practice by using the formalism to analyse 
the example of the Atwood machine. 
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example 7.8 Atwood’s Machine Using a Lagrange Multiplier 

Analyze the Atwood machine of Figure 7.6 (shown again here as Figure 7.11) 
by the method of Lagrange multipliers and using the coordinates x and y of the 
two masses. j 

In terms of the given coordinates, the Lagrangian is 

H — T — U — \rnix 1 + \m 2 y 1 + m 2 gy (7.123) j 

and the constraint equation is | 

/(x, y) = x + y = const. (7.124) 

The modified Lagrange equation (7.118) for x reads j 

9L 9/ d 9L , .. 1 

-h A— =-or m ] g-\-k = m ] x (7.125) 

dx dx dt 9x 1 1 

and that for y is 

9£ 3/ d 9£ , .. 1 

—- + X—= —— or m 2 g + k = m 2 y. (7.126) j 

9y dy dt ay 

These two equations, together with the constraint equation (7.124), are easily 
solved for the unknowns x(t), y(t ), and X(t). From (7.124) we see that y = —x, 
and then subtracting (7.126) from (7.125) we can eliminate A and arrive at the 
same result as before, j 

x = (m x - m 2 )g/(m l + m 2 ). 

I 

To better understand the two modified Lagrange equations (7.125) and 1 
(7.126), it is helpful to compare them with the two equations of the Newtonian 
solution. Newton’s second law for m x is 

m x g — F t = m x x | 



Figure 7.11 The Atwood machine again. 
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where F t is the tension in the string, and that for m 2 is 

1 m 2§ - F t = m 2y- 

| 

These are precisely the two Lagrange equations (7.125) and (7.126), with the 
Lagrange multiplier identified as the constraint force 

| 

i a. = -/•;. 

| 

[Two small comments: The minus sign occurs because both coordinates x and y 
were measured downward, whereas both tension forces are upward. In general, 
according to (7.122) the constraint force is Xdf/dx, but in this simple case, 

I df/dx = L] 

You can find some more examples of the use of Lagrange multipliers in Problems 
7.50 through 7.52. 


Principal Definitions and Equations of Chapter 7 _ 

The Lagrangian 

The Lagrangian L of a conservative system is defined as 

L = T - [/, [Eq. (7.3)] 

where T and U are respectively the kinetic and potential energies. 


Generalized Coordinates 

The n parameters q i9 • • •, q n are generalized coordinates for an N-particle system if 
every particle’s position r a can be expressed as a function of q h • • •, q n (and possibly 
the time t) and vice versa, and if n is the smallest number that allows the system to 
be described in this way. [Eqs. (7.34) & (7.35)] 

If n < 3N (in three dimensions) the system is said to be constrained. The coor¬ 
dinates q h • • •, q n are said to be natural if the functional relationships of the r a to 
q h • • •, q n are independent of time. The number of degrees of freedom of a system 
is the number of coordinates that can be independently varied. If the number of de¬ 
grees of freedom is equal to the number of generalized coordinates (in some sense the 
“normal” state of affairs), the system is said to be holonomic. [Section 7.3] 


Lagrange’s Equations 

For any holonomic system, Newton’s second law is equivalent to the n Lagrange 
equations 

9<C d 9L 
< dq [ dt d(fa 


[/ = 1, • • •, n\ [Sections 7.3 & 7.4] 
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and the Lagrange equations are in turn equivalent to Hamilton’s principle — a fact we 
used only to prove the Lagrange equations. [Eq. (7.8)] 


Generalized Momenta and Ignorable Coordinates 


The ith generalized momentum p { is defined to be the derivative 



If dh/dqi = 0, then we say the coordinate q t is ignorable and the corresponding 
generalized momentum is constant. [Section 7.6] 


The Hamiltonian 

The Hamiltonian is defined as 

n 

^ = [Eq. (7.90)] 

1 = 1 

If 3£/3 1 = 0, then TC is conserved; if the coordinates q h • • •, q n are natural, is just 
the energy of the system. 


Lagrangian for a Charge in an Electromagnetic Field 

The Lagrangian for a charge q in an electromagnetic field is 

£(r, r, t) = \mr 2 - q(V - r-A). [Eq. (7.103)] 


Problems for Chapter 7 _ 

Stars indicate the approximate level of difficulty, from easiest (*) to most difficult (★★★). 


section 7.1 Lagrange’s Equations for Unconstrained Motion 

7.1 ★ Write down the Lagrangian for a projectile (subject to no air resistance) in terms of its Cartesian 
coordinates (x, y, z), with z measured vertically upward. Find the three Lagrange equations and show 
that they are exactly what you would expect for the equations of motion. 

7.2 ★ Write down the Lagrangian for a one-dimensional particle moving along the * axis and subject 
to a force F = —kx (with k positive). Find the Lagrange equation of motion and solve it. 

7.3 ★ Consider a mass m moving in two dimensions with potential energy U(x,y) = \kr 2 , where 
r 2 = x 2 + y 2 . Write down the Lagrangian, using coordinates x and y 9 and find the two Lagrange 
equations of motion. Describe their solutions. [This is the potential energy of an ion in an “ion trap,” 
which can be used to study the properties of individual atomic ions.] 
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7.4 ★ Consider a mass m moving in a frictionless plane that slopes at an angle a with the horizontal. 
Write down the Lagrangian in terms of coordinates x, measured horizontally across the slope, and y, 
measured down the slope. (Treat the system as two-dimensional, but include the gravitational potential 
energy.) Find the two Lagrange equations and show that they are what you should have expected. 

7.5 ★ Find the components of V/(r, 0) in two-dimensional polar coordinates. [Hint: Remember that 
the change in the scalar / as a result of an infinitesimal displacement dr is df = V/ • dr.] 

7.6 ★ Consider two particles moving unconstrained in three dimensions, with potential energy 
U(r h r 2 ). (a) Write down the six equations of motion obtained by applying Newton’s second law 
to each particle, (b) Write down the Lagrangian £(r 1? r 2 , r 1? r 2 ) = T — U and show that the six La¬ 
grange equations are the same as the six Newtonian equations of part (a). This establishes the validity of 
Lagrange’s equations in rectangular coordinates, which in turn establishes Hamilton’s principle. Since 
the latter is independent of coordinates, this proves Lagrange’s equations in any coordinate system. 

7.7 ★ Do Problem 7.6, but for N particles moving unconstrained in three dimensions (in which case 
there are 3N equations of motion). 

7.8 ★★ (a) Write down the Lagrangian £ (x 1? x 2 , x h x 2 ) for two particles of equal masses, m l = m 2 = m, 
confined to the x axis and connected by a spring with potential energy U = {kx 2 . [Here x is the 
extension of the spring, x = (x x — x 2 — /), where / is the spring’s unstretched length, and I assume 
that mass 1 remains to the right of mass 2 at all times.] (b) Rewrite £ in terms of the new variables 
X = t(x L -f- x 2 ) (the CM position) and x (the extension), and write down the two Lagrange equations 
for X and x. (c) Solve for X(t) and x(0 and describe the motion. 

section 7.3 Constrained Systems in General 

7.9 ★ Consider a bead that is threaded on a rigid circular hoop of radius R lying in the xy plane with its 
center at O , and use the angle 0 of two-dimensional polar coordinates as the one generalized coordinate 
to describe the bead’s position. Write down the equations that give the Cartesian coordinates (x, y) in 
terms of 0 and the equation that gives the generalized coordinate 0 in terms of (x, y). 

7.10 ★ A particle is confined to move on the surface of a circular cone with its axis on the z axis, vertex at 
the origin (pointing down), and half-angle a . The particle’s position can be specified by two generalized 
coordinates, which you can choose to be the coordinates (p, 0) of cylindrical polar coordinates. Write 
down the equations that give the three Cartesian coordinates of the particle in terms of the generalized 
coordinates (p, 0) and vice versa. 

7.11 ★ Consider the pendulum of Figure 7.4, suspended inside a railroad car, but suppose that the car 
is oscillating back and forth, so that the point of suspension has position x s = A cos cot , y s = 0. Use the 
angle 0 as the generalized coordinate and write down the equations that give the Cartesian coordinates 
of the bob in terms of 0 and vice versa. 

section 7.4 Proof of Lagrange’s Equations with Constraints 

7.12 ★ Lagrange’s equations in the form discussed in this chapter hold only if the forces (at least the 
nonconstraint forces) are derivable from a potential energy. To get an idea how they can be modified 
to include forces like friction, consider the following: A single particle in one dimension is subject to 
various conservative forces (net conservative force = F = —dU/dx) and a nonconservative force (let’s 
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Figure 7.12 Problem 7.14 


call it F fric ). Define the Lagrangian as £ = T — U and show that the appropriate modification is 


ac 

dx 


+ ^fric — 


d dL 


dt dx 


7.13 ★★ In Section 7.4 [Equations (7.41) through (7.51)], I proved Lagrange’s equations for a single 
particle constrained to move on a two-dimensional surface. Go through the same steps to prove 
Lagrange’s equations for a system consisting of two particles subject to various unspecified constraints. 
[Hint: The net force on particle 1 is the sum of the total constraint force F^ str and the total nonconstraint 
force F 1? and likewise for particle 2. The constraint forces come in many guises (the normal force of a 
surface, the tension force of a string tied between the particles, etc.), but it is always true that the net 
work done by all constraint forces in any displacement consistent with the constraints is zero — this is 
the defining property of constraint forces. Meanwhile, we take for granted that the nonconstraint forces 
are derivable from a potential energy U(r h r 2 , t); that is, = — Vj£/ and likewise for particle 2. Write 
down the difference 8S between the action integral for the right path given by and r 2 (t) and any 
nearby wrong path given by r { (t) + € 1 (t) and r 2 (t) + e 2 (t). Paralleling the steps of Section 7.4, you 
can show that SS is given by an integral analogous to (7.49), and this is zero by the defining property 
of constraint forces.] 


section 7.5 Examples of Lagrange’s Equations 

7.14 ★ Figure 7.12 shows a crude model of a yoyo. A massless string is suspended vertically from 
a fixed point and the other end is wrapped several times around a uniform cylinder of mass m and 
radius R. When the cylinder is released it moves vertically down, rotating as the string unwinds. 
Write down the Lagrangian, using the distance x as your generalized coordinate. Find the Lagrange 
equation of motion and show that the cylinder accelerates downward with x = 2g/3. [Hints: You need 
to remember from your introductory physics course that the total kinetic energy of a body like the yoyo 
is T = i m v 2 + \Ico 2 , where v is the velocity of the center of mass, / is the moment of inertia (for a 
unform cylinder, / = \mR 2 ) and co is the angular velocity about the CM. You can express co in terms 
of x.] 

7.15 ★ A mass mi rests on a frictionless horizontal table and is attached to a massless string. The string 
runs horizontally to the edge of the table, where it passes over a massless, frictionless pulley and then 
hangs vertically down. A second mass m 2 is now attached to the bottom end of the string. Write down 
the Lagrangian for the system. Find the Lagrange equation of motion, and solve it for the acceleration of 
the blocks. For your generalized coordinate, use the distance x of the second mass below the tabletop. 
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7.16 ★ Write down the Lagrangian for a cylinder (mass m, radius R , and moment of inertia /) that 
rolls without slipping straight down an inclined plane which is at an angle a from the horizontal. Use 
as your generalized coordinate the cylinder’s distance x measured down the plane from its starting 
point. Write down the Lagrange equation and solve it for the cylinder’s acceleration x. Remember that 
T = \mv 2 + \lco 2 , where v is the velocity of the center of mass and co is the angular velocity. 

7.17 ★ Use the Lagrangian method to find the acceleration of the Atwood machine of Example 7.3 
(page 255) including the effect of the pulley’s having moment of inertia /. (The kinetic energy of the 
pulley is \Ico 2 , where go is its angular velocity.) 

7.18 ★ A mass m is suspended from a massless string, the other end of which is wrapped several times 
around a horizontal cylinder of radius R and moment of inertia /, which is free to rotate about a fixed 
horizontal axle. Using a suitable coordinate, set up the Lagrangian and the Lagrange equation of motion, 
and find the acceleration of the mass m. [The kinetic energy of the rotating cylinder is ±Ico 2 .] 

7.19 ★ In Example 7.5 (page 258) the two accelerations are given by Equations (7.66) and (7.67). 
Check that the acceleration of the block is given correctly in the limit M -» 0. [You need to find the 
components of this acceleration relative to the table .] 

7.20 ★ A smooth wire is bent into the shape of a helix, with cylindrical polar coordinates p = R and 
z = Xcjo, where R and X are constants and the z axis is vertically up (and gravity vertically down). Using 
z as your generalized coordinate, write down the Lagrangian for a bead of mass m threaded on the wire. 
Find the Lagrange equation and hence the bead’s vertical acceleration z. In the limit that R —> 0, what 
is z? Does this make sense? 

7.21 ★ The center of a long frictionless rod is pivoted at the origin, and the rod is forced to rotate in 
a horizontal plane with constant angular velocity co. Write down the Lagrangian for a bead threaded 
on the rod, using r as your generalized coordinate, where r, </> are the polar coordinates of the bead. 
(Notice that 0 is not an independent variable since it is fixed by the rotation of the rod to be 0 = cot.) 
Solve Lagrange’s equation for r(t). What happens if the bead is initially at rest at the origin? If it is 
released from any point r 0 > 0, show that r(t) eventually grows exponentially. Explain your results in 
terms of the centrifugal force mco 2 r. 

7.22 ★ Using the usual angle 0 as generalized coordinate, write down the Lagrangian for a simple 
pendulum of length / suspended from the ceiling of an elevator that is accelerating upward with 
constant acceleration a. (Be careful when writing 7; it is probably safest to write the bob’s velocity 
in component form.) Find the Lagrange equation of motion and show that it is the same as that for a 
normal, nonaccelerating pendulum, except that g has been replaced by g -f a. In particular, the angular 
frequency of small oscillations is ^J(g + a)/1 . 

7.23 ★ A small cart (mass m) is mounted on rails inside a large cart. The two are attached by a spring 
(force constant k) in such a way that the small cart is in equilibrium at the midpoint of the large. The 
distance of the small cart from its equilibrium is denoted x and that of the large one from a fixed 
point on the ground is X, as shown in Figure 7.13. The large cart is now forced to oscillate such that 
X = A coscot, with both A and co fixed. Set up the Lagrangian for the motion of the small cart and 
show that the Lagrange equation has the form 

x + co 2 x = B cos cot 

where co 0 is the natural frequency co Q = V k/m and B is a constant. This is the form assumed in Section 
5.5, Equation (5.57), for driven oscillations (except that we are here ignoring damping). Thus the system 
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described here would be one way to realize the motion discussed there. (We could fill the large cart 
with molasses to provide some damping.) 

7.24 ★ We saw in Example 7.3 (page 255) that the acceleration of the Atwood machine is x = 
(m { — m 2 )g/(m l + m 2 ). It is sometimes claimed that this result is “obvious” because, it is said, the 
effective force on the system is (m x — m 2 )g and the effective mass is (raj + m 2 ). This is not, perhaps, 
all that obvious, but it does emerge very naturally in the Lagrangian approach. Recall that Lagrange’s 
equation can be thought of as [Equation (7.17)] 

(generalized force) = (rate of change of generalized momentum). 

Show that for the Atwood machine the generalized force is (m l — m 2 )g and the generalized momentum 
(m x + m 2 )x. Comment. 

7.25 ★ Prove that the potential energy of a central force F = —kr n r (with n — 1) is U = kr n+l /(n + 
1). In particular, if n = 1, then F = —kr and U = \kr 2 . 

7.26 ★ In Example 7.7 (page 264), we saw that the bead on a spinning hoop can make small oscillations 
about any of its stable equilibrium points. Verify that the oscillation frequency Q' defined in (7.79) is 
equal to j co 2 — ( g/coR ) 2 as claimed in (7.80). 

7.27 ★★ Consider a double Atwood machine constructed as follows: A mass 4m is suspended from a 
string that passes over a massless pulley on frictionless bearings. The other end of this string supports 
a second similar pulley, over which passes a second string supporting a mass of 3m at one end and m 
at the other. Using two suitable generalized coordinates, set up the Lagrangian and use the Lagrange 
equations to find the acceleration of the mass 4m when the system is released. Explain why the top 
pulley rotates even though it carries equal weights on each side. 

7.28 ★★ A couple of points need checking from Example 7.6 (page 260). (a) From the point of view of 
a noninertial frame rotating with the hoop, the bead is subject to the force of gravity and a centrifugal 
force mco 2 p (in addition to the constraint force, which is the normal force of the wire). Verify that at the 
equilibrium points given by (7.71), the tangential components of these two forces balance one another. 
(A free-body diagram will help.) (b) Verify that the equilibrium point at the top (0 = n) is unstable, 
(c) Verify that the equilibrium at the second point given by (7.71) (the one on the left, with 0 negative) 
is stable. 

7.29 ★★ Figure 7.14 shows a simple pendulum (mass m, length /) whose point of support P is attached 
to the edge of a wheel (center 0, radius R ) that is forced to rotate at a fixed angular velocity co. At 
t = 0, the point P is level with O on the right. Write down the Lagrangian and find the equation of 
motion for the angle 0. [Hint: Be careful writing down the kinetic energy T. A safe way to get the 
velocity right is to write down the position of the bob at time t, and then differentiate.] Check that your 
answer makes sense in the special case that co — 0. 
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Figure 7.14 Problem 7.29 


7.30 ★★ Consider the pendulum of Figure 7.4, suspended inside a railroad car that is being forced to 
accelerate with a constant acceleration a. (a) Write down the Lagrangian for the system and the equation 
of motion for the angle 0. Use a trick similar to the one used in Equation (5.11) to write the combination 
of sin 0 and cos 0 as a multiple of sin (0 4-/3). (b) Find the equilibrium angle 0 at which the pendulum 
can remain fixed (relative to the car) as the car accelerates. Use the equation of motion to show that 
this equilbrium is stable. What is the frequency of small oscillations about this equilibrium position? 
(We shall find a much slicker way to solve this problem in Chapter 9, but the Lagrangian method does 
give a straightforward route to the answer.) 

7.31 ★★ A simple pendulum (mass M and length L) is suspended from a cart (mass m) that can oscillate 
on the end of a spring of force constant k, as shown in Figure 7.15. (a) Write the Lagrangian in terms 
of the two generalized coordinates x and 0, where x is the extension of the spring from its equilibrium 
length. (Read the hint in Problem 7.29.) Find the two Lagrange equations. (Warning: They’re pretty 
ugly!) (b) Simplify the equations to the case that both x and 0 are small. (They’re still pretty ugly, and 
note, in particular, that they are still coupled ; that is, each equation involves both variables. Nonetheless, 
we shall see how to solve these equations in Chapter 11 — see particularly Problem 11.19.) 



Figure 7.15 Problem 7.31 


7.32 ★★ Consider the cube balanced on a cylinder as described in Example 4.7 (page 130). Assuming 
that b < r, use the Lagrangian approach to find the angular frequency of small oscillations about the 
top. The simplest procedure is to make the small-angle approximations to £ before you differentiate to 
get Lagrange’s equation. As usual, be careful in writing down the kinetic energy; this is \ (mv 2 4- 10 2 ), 
where v is the speed of the CM and I is the moment of inertia about the CM (2mb 2 /3). The safe way 
to find v is to write down the coordinates of the CM and then differentiate. 
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7.33 ★★ A bar of soap (mass m) is at rest on a frictionless rectangular plate that rests on a horizontal 
table. At time t = 0, I start raising one edge of the plate so that the plate pivots about the opposite 
edge with constant angular velocity co , and the soap starts to slide toward the downhill edge. Show that 
the equation of motion for the soap has the form x — co 2 x = —g sin cot , where x is the soap’s distance 
from the downhill edge. Solve this for x(t), given that x{0) — x 0 . [You’ll need to use the method used 
to solve Equation (5.48). You can easily solve the homogeneous equation; for a particular solution try 
x = A sin cot and solve for A.] 

7.34 ★★ Consider the well-known problem of a cart of mass m moving along the x axis attached to a 
spring (force constant k), whose other end is held fixed (Figure 5.2). If we ignore the mass of the spring 
(as we almost always do) then we know that the cart executes simple harmonic motion with angular 
frequency co = +Jk/m . Using the Lagrangian approach, you can find the effect of the spring’s mass 
M, as follows: (a) Assuming that the spring is uniform and stretches uniformly, show that its kinetic 
energy is | Mx 2 . (As usual x is the extension of the spring from its equilibrium length.) Write down 
the Lagrangian for the system of cart plus spring. {Note: The potential energy is still \kx 2 .) (b) Write 
down the Lagrange equation and show that the cart still executes SHM but with angular frequency 
co — yj.k/{m 4- M/3); that is, the effect of the spring’s mass M is just to add Mj 3 to the mass of the 
cart. 



Figure 7.16 Problem 7.35 


7.35 ★★ Figure 7.16 is a bird’s-eye view of a smooth horizontal wire hoop that is forced to rotate at a 
fixed angular velocity co about a vertical axis through the point A. A bead of mass m is threaded on the 
hoop and is free to move around it, with its position specified by the angle 0 that it makes at the center 
with the diameter AB. Find the Lagrangian for this system using (j> as your generalized coordinate. 
(Read the hint in Problem 7.29.) Use the Lagrange equation of motion to show that the bead oscillates 
about the point B exactly like a simple pendulum. What is the frequency of these oscillations if their 
amplitude is small? 

7.36 ★★★ A pendulum is made from a massless spring (force constant k and unstretched length l 0 ) that 
is suspended at one end from a fixed pivot O and has a mass m attached to its other end. The spring 
can stretch and compress but cannot bend, and the whole system is confined to a single vertical plane, 
(a) Write down the Lagrangian for the pendulum, using as generalized coordinates the usual angle <j> 
and the length r of the spring, (b) Find the two Lagrange equations of the system and interpret them 
in terms of Newton’s second law, as given in Equation (1.48). (c) The equations of part (b) cannot be 
solved analytically in general. However, they can be solved for small oscillations. Do this and describe 
the motion. [Hint: Let / denote the equilibrium length of the spring with the mass hanging from it and 
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write r = / + e. “Small oscillations” involve only small values of e and 0, so you can use the small- 
angle approximations and drop from your equations all terms that involve powers of 6 or 0 (or their 
derivatives) higher than the first power (also products of e and 0 or their derivatives). This dramatically 
simplifies and uncouples the equations.] 

7.37 ★★★ Two equal masses, mj = m 2 — ra, are joined by a massless string of length L that passes 
through ^ hole in a frictionless horizontal table. The first mass slides on the table while the second 
hangs below the table and moves up and down in a vertical line, (a) Assuming the string remains taut, 
write down the Lagrangian for the system in terms of the polar coordinates (r, 0) of the mass on the table, 
(b) Find the two Lagrange equations and interpret the 0 equation in terms of the angular momentum t 
of the first mass, (c) Express 0 in terms of i and eliminate 0 from the r equation. Now use the r equation 
to find the value r = r 0 at which the first mass can move in a circular path. Interpret your answer in 
Newtonian terms, (d) Suppose the first mass is moving in this circular path and is given a small radial 
nudge. Write r(t ) = r 0 + e{t) and rewrite the r equation in terms of e(t) dropping all powers of e(t) 
higher than linear. Show that the circular path is stable and that r(t) oscillates sinusoidally about r 0 . 
What is the frequency of its oscillations? 

7.38 ★★★ A particle is confined to move on the surface of a circular cone with its axis on the vertical 
z axis, vertex at the origin (pointing down), and half-angle a. (a) Write down the Lagrangian L in 
terms of the spherical polar coordinates r and 0. (b) Find the two equations of motion. Interpret the 0 
equation in terms of the angular momentum i z , and use it to eliminate 0 from the r equation in favor 
of the constant i z . Does your r equation make sense in the case that l z — 0? Find the value r 0 of r 
at which the particle can remain in a horizontal circular path, (c) Suppose that the particle is given a 
small radial kick, so that r(t) = r Q -f e(t), where e(t) is small. Use the r equation to decide whether 
the circular path is stable. If so, with what frequency does r oscillate about r G ? 

7.39 ★★★ (a) Write down the Lagrangian for a particle moving in three dimensions under the influence 
of a conservative central force with potential energy U (r), using spherical polar coordinates (r, 0,0). 
(b) Write down the three Lagrange equations and explain their significance in terms of radial accelera¬ 
tion, angular momentum, and so forth. (The 0 equation is the tricky one, since you will find it implies 
that the 0 component of i varies with time, which seems to contradict conservation of angular mo¬ 
mentum. Remember, however, that is the component of i in a variable direction.) (c) Suppose that 
initially the motion is in the equatorial plane (that is, 6 0 = tt/2 and 0 O = 0). Describe the subsequent 
motion, (d) Suppose instead that the initial motion is along a line of longitude (that is, 0 O = 0). Describe 
the subsequent motion. 

7.40 ★★★ The “spherical pendulum” is just a simple pendulum that is free to move in any sideways 
direction. (By contrast a “simple pendulum” — unqualified — is confined to a single vertical plane.) 
The bob of a spherical pendulum moves on a sphere, centered on the point of support with radius r = R, 
the length of the pendulum. A convenient choice of coordinates is spherical polars, r, 0, 0, with the 
origin at the point of support and the polar axis pointing straight down. The two variables 0 and 0 make 
a good choice of generalized coordinates, (a) Find the Lagrangian and the two Lagrange equations, 
(b) Explain what the 0 equation tells us about the z component of angular momentum t z . (c) For the 
special case that 0 = const, describe what the 0 equation tells us. (d) Use the 0 equation to replace 0 
by t z in the 0 equation and discuss the existence of an angle 0 O at which 0 can remain constant. Why 
is this motion called a conical pendulum? (e) Show that if 0 = 0 O + e, with e small, then 0 oscillates 
about 0 O in harmonic motion. Describe the motion of the pendulum’s bob. 

7.41 ★★★ Consider a bead of mass m sliding without friction on a wire that is bent in the shape of a 
parabola and is being spun with constant angular velocity co about its vertical axis, as shown in Figure 
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7.17. Use cylindrical polar coordinates and let the equation of the parabola be z = kp 2 . Write down the 
Lagrangian in terms of p as the generalized coordinate. Find the equation of motion of the bead and 
determine whether there are positions of equilibrium, that is, values of p at which the bead can remain 
fixed, without sliding up or down the spinning wire. Discuss the stability of any equilibrium positions 
you find. 



Figure 7.17 Problem 7.41 


7.42 ★★★ [Computer] In Example 7.7 (page 264), we saw that the bead on a spinning hoop can make 
small oscillations about its nonzero stable equilibrium points that are approximately sinusoidal, with 
frequency Q' = y/co 2 — ( g/coR ) 2 as in (7.80). Investigate how good this approximation is by solving 
the equation of motion (7.73) numerically and then plotting both your numerical solution and the 
approximate solution 6(t) = 0 o -F Acos(Q't — 8) on the same graph. Use the following numbers: 
g = R = 1 and co 2 = 2, and initial conditions 0(0) = 0 and 0(0) = 0 O + € 0 , where e 0 = 1°. Repeat 
with € 0 = 10°. Comment on your results. 

7.43 ★★★ [Computer] Consider a massless wheel of radius R mounted on a frictionless horizontal axis. 
A point mass M is glued to the edge, and a massless string is wrapped several times around the perimeter 
and hangs vertically down with a mass m suspended from its bottom end. (See Figure 4.28.) Initially I 
am holding the wheel with M vertically below the axle. At t = 0,1 release the wheel, and m starts to fall 
vertically down, (a) Write down the Lagrangian £ = T — U as a function of the angle 0 through which 
the wheel has turned. Find the equation of motion and show that, provided m < M, there is one position 
of stable equilibrium, (b) Assuming m < M, sketch the potential energy U (0) for — ir < 0 < 47r and 
use your graph to explain the equilibrium position you found, (c) Because the equation of motion 
cannot be solved in terms of elementary functions, you are going to solve it numerically. This requires 
that you choose numerical values for the various parameters. Take M = g = R = 1 (this amounts to 
a convenient choice of units) and m =0.7. Before solving the equation make a careful plot of U (0) 
against 0 and predict the kind of motion expected when M is released from rest at 0 = 0. Now solve 
the equation of motion for 0 < t < 20 and verify your prediction, (d) Repeat part (c), but with m = 0.8. 

7.44 ★★★ [Computer] If you haven’t already done so, do Problem 7.29. One might expect that the 
rotation of the wheel would have little effect on the pendulum, provided the wheel is small and rotates 
slowly, (a) Verify this expectation by solving the equation of motion numerically, with the following 
numbers: Take g and l to be 1. (This means that the natural frequency y/g/1 of the pendulum is also 
1.) Take co = 0.1, so that the wheel’s rotational frequency is small compared to the natural frequency 
of the pendulum; and take the radius R = 0.2, significantly less than the length of the pendulum. As 
initial conditions take 0 = 0.2 and 0 = 0 at t = 0, and make a plot of your solution 0 ( t) for 0 < t < 20. 
Your graph should look very like the sinusoidal oscillations of an ordinary simple pendulum. Does the 
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period look correct? (b) Now plot 0(0 for 0 < t < 100 and notice that the rotating support does make 
a small difference, causing the amplitude of the oscillations to grow and shrink periodically. Comment 
on the period of these small fluctuations. 


section 7.8 More About Conservation Laws* 


7.45 ★★ (a) Verify that the coefficients A /; - in the important expression (7.94) for the kinetic energy of 
any “natural” system are symmetric; that is, A t j = Aj t . (b) Prove that for any n variables v 1? • • •, v n 

-^Y, A Jk v J v k= 2 Y2 A u v J- 

1 M j 

[Hint: Start with the case that n = 2, for which you can write out the sums in full. Notice that you 
need the result of part (a).] This identity is useful in many areas of physics; we needed it to prove the 
expression (7.96) for the generalized momentum p { . 

7.46 ★★ Noether’s theorem asserts a connection between invariance principles and conservation laws. 
In Section 7.8 we saw that translational invariance of the Lagrangian implies conservation of total linear 
momentum. Here you will prove that rotational invariance of L implies conservation of total angular 
momentum. Suppose that the Lagrangian of an A-particle system is unchanged by rotations about a 
certain symmetry axis, (a) Without loss of generality, take this axis to be the z axis, and show that 
the Lagrangian is unchanged when all of the particles are simultaneously moved from (r a , 0 a , 0 a ) to 
(r a , 0 a , (p a + e) (same 6 for all particles). Hence show that 


£ 


dL 

d(j) a 


= 0 . 


(b) Use Lagrange’s equations to show that this implies that the total angular momentum L z about the 
symmetry axis is constant. In particular, if the Lagrangian is invariant under rotations about all axes, 
then all components of L are conserved. 

7.47 ★★★ In Chapter 4 (at the end of Section 4.7) I claimed that, for a system with one degree of freedom, 
positions of stable equilibrium “normally” correspond to minima of the potential energy U ( q ). Using 
Lagrangian mechanics, you can now prove this claim, (a) Consider a one-degree system of N particles 
with positions r a = r a (q), where q is the one generalized coordinate and the transformation between 
r and q does not depend on time; that is, q is what we have now agreed to call “natural.” (This is the 
meaning of the qualification “normally” in the statement of the claim. If the transformation depends 
on time, then the claim is not necessarily true.) Prove that the KE has the form T = \Aq 2 , where 
A = A{q) > 0 may depend on q but not on q. [This corresponds exactly to the result (7.94) for n 
degrees of freedom. If you have trouble with the proof here, review the proof there.] Show that the 
Lagrange equation of motion has the form 


A(q)q = 


dU 

dq 


IdA . 2 

- - q 

2 dq 


(b) A point q 0 is an equilibrium point if, when the system is placed at q Q with q — 0, it remains there. 
Show that q Q is an equilibrium point if and only if dU I dq — 0. (c) Show that the equilibrium is stable 
if and only if U is minimum at q 0 . (d) If you did Problem 7.30, show that the pendulum of that problem 
does not satisfy the conditions of this problem and that the result proved here is false for that system. 
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section 7.9 Lagrange’s Equations for Magnetic Forces* 

7.48 ★★ Let F = F(q h • • •, q n ) be any function of the generalized coordinates (q h • • •, q n ) of a system 
with Lagrangian L(q h • • •, q n , q h • • •, q n , t). Prove that the two Lagrangians £ and £/ = £ + dF/dt 
give exactly the same equations of motion. 

7.49 ★★ Consider a particle of mass m and charge q moving in a uniform constant magnetic field B in 
the z direction, (a) Prove that B can be written as B = V x A with A = |B x r. Prove equivalently 

/V 

that in cylindrical polar coordinates, A — \Bp$. (b) Write the Lagrangian (7.103) in cylindrical polar 
coordinates and find the three corresponding Lagrange equations, (c) Describe in detail those solutions 
of the Lagrange equations in which p is a constant. 

section 7 .io Lagrange Multipliers and Constraint Forces* 

7.50 ★ A mass m j rests on a frictionless horizontal table. Attached to it is a string which runs hori¬ 
zontally to the edge of the table, where it passes over a frictionless, small pulley and down to where 
it supports a mass m 2 . Use as coordinates x and y the distances of m l and m 2 from the pulley. These 
satisfy the constraint equation /(x, y) = x + y = const. Write down the two modified Lagrange equa¬ 
tions and solve them (together with the constraint equation) for x, y, and the Lagrange multiplier A. 
Use (7.122) (and the corresponding equation in y) to find the tension forces on the two masses. Verify 
your answers by solving the problem by the elementary Newtonian approach. 

7.51 ★ Write down the Lagrangian for the simple pendulum of Figure 7.2 in terms of the rectangular 
coordinates x and y. These coordinates are constrained to satisfy the constraint equation /(x, y) = 
^jx 2 + y 2 — l. (a) Write down the two modified Lagrange equations (7.118) and (7.119). Comparing 
these with the two components of Newton’s second law, show that the Lagrange multiplier is (minus) 
the tension in the rod. Verify Equation (7.122) and the corresponding equation in y. (b) The constraint 
equation can be written in many different ways. For example we could have written /'(x, y) = 
x 2 -p y 2 = l 2 . Check that using this function would have given the same physical results. 

7.52 ★ The method of Lagrange multipliers works perfectly well with non-Cartesian coordinates. 
Consider a mass m that hangs from a string, the other end of which is wound several times around 
a wheel (radius R , moment of inertia I) mounted on a frictionless horizontal axle. Use as coordinates 
for the mass and the wheel x, the distance fallen by the mass, and 0, the angle through which the 
wheel has turned (both measured from some convenient reference position). Write down the modified 
Lagrange equations for these two variables and solve them (together with the constraint equation) for 
x and </> and the Lagrange multiplier. Write down Newton’s second law for the mass and wheel, and 
use them to check your answers for x and <j>. Show that Xdf/dx is indeed the tension force on the mass. 
Comment on the quantity A 9// 3<j>. 



